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Abstract 


The  original  version  of  the  4SIGHT  computer  program  was  written  to  facilitate  the 
performance  assessment  of  buried  concrete  structures  based  deterministic  calculations 
of  service  life,  using  data  supplied  by  the  user.  In  an  effort  to  ensure  reliability  in  the 
results,  a means  of  validation  was  required.  Since  the  results  from  the  calculations  depend 
strongly  on  the  diffusive  and  advective  transport  of  ionic  species,  validating  the  transport 
portion  of  4SIGHT  was  paramount. 

The  validation  of  the  transport  portion  of  4SIGHT  used  both  reference  and  laboratory 
data.  The  reference  data  include  electrolyte  densities,  activity  coefficients,  and  binary 
salt  diffusion  coefficients.  The  laboratory  experiments  included  diffusion  experiments  per- 
formed on  relatively  inert  ceramic  frits,  using  a range  of  ionic  species  and  concentrations. 
The  laboratory  experiments  were  also  used  to  confirm  the  fundamental  assumption  that 
multi-species  transport  in  non-reacting  systems  can  be  characterized  by  two  numbers, 
the  porosity  and  the  formation  factor,  regardless  of  the  number  of  diffusing  species. 

A default  analysis  for  estimating  the  spacing  and  thickness  of  flexural  and  drying 
shrinkage  cracks  was  added  to  4SIGHT.  The  flexural  cracks  are  calculated  from  a rudi- 
mentary structural  analysis  that  assumes  the  roof  of  the  vault  acts  as  a one-way  slab. 
The  drying  shrinkage  cracks  are  estimated,  in  large  part,  from  the  total  water  content 
of  the  concrete  mixture,  the  exposure  environment,  and  the  geometry  of  the  slab.  The 
subsequent  predicted  crack  widths,  crack  spacings,  and  location  of  the  neutral  axis,  are 
used  to  predict  the  bulk  permeability  of  the  concrete  structure  as  a function  of  depth; 
the  portion  above  the  neutral  axis  being  uncracked,  the  portion  below  the  neutral  axis 
being  cracked. 

Another  addition  to  4SIGHT  was  the  means  to  perform  a Monte  Carlo  calculation, 
based  on  parameter  uncertainty,  of  the  duration  of  useful  service.  The  service  life  calcu- 
lation is  repeated  multiple  times  using  a random  number  generator  to  calculate  random 
deviates  for  the  physical  parameters.  The  results  can  then  be  used  to  make  a probabilistic 
statement  about  the  performance  of  the  concrete  structure. 

The  resulting  computer  program  is  state-of-the-art  attempt  to  characterize  both  the 
response  of  a concrete  to  its  environment.  However,  there  can  be  no  assurances  that  the 
model  is  exact  for  every  scenario  possible,  regardless  of  parameter  uncertainty.  Therefore, 
in  cases  that  concern  public  safety,  the  only  way  to  confirm  the  response  of  the  concrete 
to  its  environment  is  through  continuous  monitoring  and  periodic  testing.  The  computer 
program,  in  these  circumstances,  can  also  serve  as  a useful  tool  in  establishing  a testing 
monitoring  protocol.  The  spatial  and  temporal  distribution  of  the  monitoring  and  testing 
protocol  can  be  established  using  the  computer  program  as  a guide. 
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1 Introduction 


The  4SIGHT  computer  program,  first  published  in  1995  [1],  was  developed  as  a resource 
for  estimating  the  service  life  of  new  underground  concrete  structures.  The  program 
was  the  first  to  use  combined  numerical  models  for  ion  transport,  chemical  reaction,  and 
subsequent  changes  to  transport  coefficients  to  model  the  response  of  a concrete  structure 
to  its  environment.  Although  the  individual  components  of  4SIGHT  duplicated  previous 
work,  the  unique  ensemble  allowed  4SIGHT  to  model  the  synergistic  degradation  effects 
that  may  occur.  While  this  approach  is  physically  sound,  the  model  remained  to  be 
validated  in  order  to  ensure  confidence  in  the  results. 

The  initial  validation  concept  was  to  use  samples  from  field  concrete  and  compare  the 
prediction  of  4SIGHT  to  the  actual  response  of  the  concrete.  While  conceptually  sound, 
there  were  unforeseen  pitfalls  to  such  an  approach.  Since  4SIGHT  was  to  be  used  to 
predict  the  performance  over  decades,  and  possibly  centuries,  comparable  concrete  data 
were  needed.  While  such  concretes  exist,  the  documentation  of  these  structure  is  either 
scarce  or  non-existent.  Further,  the  environmental  conditions  that  the  concretes  were 
subjected  to  were  also  difficult  to  quantify.  These  undetermined  parameters  were  either 
estimated  using  other  data  within  the  report  or  left  as  a free  parameter,  and  adjusted  until 
the  4SIGHT  prediction  agreed  with  the  held  concrete  response.  Since  there  were  always 
a number  of  input  parameters  that  were  undetermined,  the  held  validation  procedure 
could  only  be  used  to  determine  if  4SIGHT  was  unreasonable;  the  calculation  could  not 
be  used  to  validate  the  program.  Although  the  results  of  these  trials  were  favorable,  they 
were  by  no  means  a validation  of  4SIGHT. 

The  held  performance  validation  concept,  which  has  been  used  elsewhere  [2],  is  a 
common  approach  in  service  life  models  for  concrete.  With  few  exceptions,  these  models 
are  empirical  and  attempt  to  predict  the  response  of  a concrete  under  certain  environ- 
mental conditions;  they  do  not  attempt  to  capture  every  detail  of  transport  and  reaction. 
Some  of  these  models  are  based  on  a transport  equation  such  as  Fick's  law.  Others  are 
phenomenological  modihcations  to  Fick’s  law  that  neglect  the  details  of  transport  and 
reaction  and  seek  to  simply  mimic  the  macroscopic  response  of  the  concrete.  Generally, 
the  model  parameters  are  estimated  from  held  data,  and  then  the  model  is  applied  to 
other  concretes. 

By  contrast,  the  approach  of  4SIGHT  is  to  capture  every  detail  of  ion  transport  and 
chemical  reaction.  In  a sense,  the  objective  is  to  simulate  reality  in  an  attempt  to  predict 
future  performance.  To  achieve  this,  the  hrst  step  is  to  fully  characterize  ion  transport 
and  reaction  within  the  concrete.  The  second  is  to  fully  characterize  the  degradation 
process,  which  includes  both  the  effects  of  reaction  on  transport  properties  and  the 
extent  of  reaction  at  which  sufficient  pressures  are  produced  so  as  to  create  mechanical 
failure  within  the  concrete.  This  approach  is  virtually  unique  among  concrete  service  life 
computer  models. 

In  order  to  validate  the  transport  portion  of  the  4SIGHT  computer  program  an  inert 
ceramic  frit  was  used  as  the  model  porous  material.  By  using  the  ceramic  frit,  the 
effects  of  transport  could  be  isolated  from  the  effects  of  reaction.  Further,  the  results  of 
impedance  spectroscopy  experiments  suggested  that  the  frits  contained  relatively  little 
surface  charge,  which  would  result  in  relatively  little  surface  binding  of  the  diffusing  ionic 
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species,  further  contributing  to  nearly  ideal  diffusive  transport  behavior.  A practical 
advantage  of  the  ceramic  over  cement  pastes  is  that  the  diffusivity  of  the  ceramic  frits 
was  relatively  large,  meaning  that  data  could  be  collected  over  a comparatively  short 
period  of  time. 

The  results  from  the  diffusion  and  the  impedance  spectroscopy  measurements  yielded 
two  important  facts:  the  diffusive  transport  of  ionic  species  within  a porous  material 
can  be  characterized  by  two  numbers,  the  porosity  and  the  formation  factor  (defined 
subsequently),  regardless  of  the  number  of  different  diffusing  species.  The  apparent 
diffusion  coefficient  that  is  often  reported  is  of  relatively  little  use,  particularly  over  long 
periods  of  time.  These  two  results  are  significant  in  that  the  vast  majority  of  concrete 
materials  transport  studies  concentrate  on  the  determination  of  an  apparent  chloride 
diffusion  coefficient  that  is  calculated  from  a simplistic  transport  equation.  The  fact  that 
formation  factor  is  far  more  useful  has  distinct  experimental  advantage:  the  formation 
factor  can  be  determined  from  an  electrical  measurement  that  can  be  performed  by 
virtually  every  concrete  materials  laboratory.  Further,  the  formation  factor  is  a concept 
that  originated  in  the  field  of  geophysical  research.  Therefore,  it  is  a technique  that 
should  also  be  familiar  to  geological  testing  laboratories. 

This  report  summarizes  the  experiments  that  were  performed  in  order  to  validate 
the  transport  portion  of  the  4SIGHT  computer  model.  Some  of  the  validation  consists  of 
comparing  the  output  from  4SIGHT  calculations  to  physical  and  chemical  handbook  data 
for  various  electrolyte  properties.  The  remaining  tests  are  laboratory  experiments  per- 
formed on  the  aforementioned  ceramic  frits.  The  details  of  the  experiments  are  published 
elsewhere,  and  the  results  are  summarized  here. 

In  addition  to  the  validation,  two  enhancements  to  4SIGHT  were  requested:  a model 
to  predict  crack  formation  within  the  concrete;  and  a means  to  perform  a Monte  Carlo 
calculation  of  the  service  life  based  on  parameter  uncertainty.  A model  for  crack  forma- 
tion was  necessary  due  to  the  dramatic  influence  of  cracks  on  advective  transport  (Darcy 
flow)  through  a porous  material.  The  Monte  Carlo  calculation  addresses  the  nature  of 
uncertainty  for  any  type  of  service  life  calculation.  No  physical  parameter  can  be  de- 
termined with  absolute  certainty,  and  in  some  cases  small  variations  in  a parameter  can 
have  a dramatic  effect  on  the  overall  performance  of  a structure.  The  Monte  Carlo  cal- 
culation was  implemented  with  a robust  random  number  generator  and  permits  the  user 
to  specify  the  parameter  uncertainty. 

The  cracking  model  considered  the  possible  causes  of  cracks  within  the  concrete  and 
divided  the  sources  of  cracking  between  “engineered”  cracks  and  construction  cracks.  The 
engineered  cracks  (flexural  and  drying  shrinkage)  can  be  reasonably  estimated  based  on 
the  structural  design,  the  materials  specification,  and  the  exposure.  By  contrast,  con- 
struction cracks  (plastic  shrinkage,  settlement,  and  thermal)  can  be  avoided  by  applying 
proper  construction  techniques;  predicting  how  well  builders  will  perform  their  duties  is 
an  open-ended  question  with  an  indeterminate  amount  of  uncertainty.  Hence,  the  crack 
prediction  was  limited  to  engineering  cracks,  and  are  estimated  from  a simple  structural 
mechanics  model  (flexural)  and  the  relationships  provided  in  the  American  Concrete 
Institute  (ACI)  Manual  of  Concrete  Practice  [3]  (drying  shrinkage). 
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2 Degradation  Validation  Examples 

The  use  of  either  literature  data  or  field  data  for  validating  4SIGHT  lacks  rigor.  In  every 
case,  a critical  piece  of  information  is  missing,  which  means  that  some  parameter  used 
in  the  4SIGHT  model  remains  as  an  independent  parameter.  This  free  parameter  can  be 
adjusted  until  the  calculated  output  from  4SIGHT  coincides  with  the  results  reported  in 
the  literature.  If  the  value  chosen  for  the  free  parameter  is  typical  for  similar  concretes, 
one  can  only  conclude  that  the  4SIGHT  model  is  not  unreasonable. 

A similar  scenario  arises  with  held  concretes.  For  a held  concrete  that  has  been  in 
service  for  decades,  many  pieces  of  data  are  missing.  In  the  rare  cases  that  the  const  ruc- 
tion was  documented,  many  pieces  of  information  required  to  make  accurate  estimates  of 
transport  properties  would  not  be  present.  Further,  for  these  concrete  structures  there 
are  virtually  no  data  of  the  environmental  conditions  the  concrete  had  been  subjected  to 
over  the  decades  of  service. 

Another  difficulty  in  the  use  of  held  concretes  was  finding  an  “olch  concrete  that 
had  been  sufficiently  documented,  yet  was  in  the  appropriate  state  of  degradation.  If  a 
concrete  structure  that  was  placed  in  the  1940s  degraded  after  two  decades  of  service, 
that  concrete  would  no  longer  exist  today.  If  the  same  concrete  was  appropriate  for  its 
intended  use,  it  would  still  be  in  service  today,  without  any  indications  of  degradation. 
Old  concrete  in  a measurable  state  of  significant  degradation  is  rare.  It  was  determined 
that  extraction  and  testing  of  held  concretes  would  not  be  cost  effective  given  the  rarity 
of  these  specimens  and  the  number  of  unresolved  parameters. 

The  following  examples  serve  as  a demonstrations  of  the  types  of  data  available  in 
the  literature.  Although  each  example  presents  all  of  the  relevant  data  presented  in  the 
report,  some  data  are  stdl  missing  and  must  be  estimated  in  some  way. 

2.1  Literature 

The  task  of  validating  4SIGHT  using  published  literature  required  hnding  suitable  reports 
of  concrete  performance  under  various  exposures.  There  are  a large  number  of  publica- 
tions characterizing  test  specimens  and  environments.  Because  4SIGHT  will  be  used  to 
predict  the  long  term  response  of  concrete,  publications  dealing  with  “old”  concrete  were 
needed.  This  fact,  alone  reduced  the  number  of  eligible  publications  dramatically. 

An  additional  complication  arises  from  the  state  of  concrete  technology  at  the  time 
the  concrete  was  produced.  Ideally,  concretes  at  least  40  years  old  would  be  desirable  for 
this  study.  The  majority  of  structures  having  this  longevity  are  typically  public  projects 
that  consumed  large  quantities  of  concrete,  the  components  of  which  varied  over  the  life 
of  the  project,  and  the  properties  of  which  were  not  characterized  thoroughly.  This  lack 
of  detailed  information  limits  present-day  characterization  of  the  concrete.  Much  of  the 
same  can  be  said  for  the  cements  used  at  the  time,  for  which  very  little  quantitative 
information  exists  regarding  the  cement  fineness,  oxide  composition,  and  mineral  com- 
position. Historical  publications  such  as  the  Building  Science  Series  [4]  published  by  the 
U.S.  Department  of  Commerce  can  only  serve  as  general  guides. 
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Figure  1:  Sample  porosity,  based  upon  pores  1 /.mi  to  10  /mi  in  size,  as  a function  of 
distance  from  the  exposure  surface.  Graph  taken  from  Tragoardh  and  Lagerblad  [5] 

2.1.1  Leaching 

In  1906,  a drinking  water  tank  was  installed  in  Uppsala,  Sweden.  The  tank  was  con- 
structed of  steel  and  lined  with  a concrete  mortar.  In  1991,  the  tank  was  demolished, 
and  the  concrete  was  analyzed  by  the  Swedish  Cement  and  Concrete  Research  Institute. 
The  results  of  the  analyses  were  published  in  an  institute  report  written  by  Tragoardh 
and  Lagerblad  [5]. 

The  water  exposure  of  the  concrete  provided  a nearly  constant  environment  to  the 
concrete.  During  its  service,  the  tank  was  periodically  drained  and  filled.  There  was 
relatively  little  evidence  of  carbonation,  suggesting  that  the  specimen  analyzed  had  been 
submerged  for  nearly  its  entire  period  of  service. 

Characterization  of  the  concrete  is  more  difficult.  No  information  was  given  regarding 
the  mixture  design.  The  only  characterization  the  report  offered  was  that  the  concen- 
tration of  calcium  hydroxide  suggested  the  water-cement  ratio  might  have  been  approx- 
imately 0.40,  which  would  have  been  quite  advanced  for  that  time. 

The  porosity  measurements  for  pores  in  the  range  1 //m  to  10  //in  have  been  repro- 
duced in  Fig.  1.  These  measurements  show  that  the  depth  of  leaching  is  no  greater  than 
7.5  mm,  and  that  the  bulk  concrete  porosity  is  3 % to  4 %.  The  corresponding  bulk  paste 
porosity  would  be  approximately  thrice  this  value,  or  9 % to  12  %.  Also,  the  measure- 
ments show  that  porosity  near  the  exposed  surface  is  nearly  three  times  greater  than  the 
bulk  porosity,  which  may  be  due,  in  part,  to  the  greater  paste  content  near  the  surface. 

Since  the  chloride  diffusion  coefficient  of  the  sample  was  not  measured,  it  had  to  be 
estimated.  The  4SIGHT  prediction  is  based  on  a simple  model  for  the  diffusion  coefficient 
(nr-s-1)  that  is  only  a function  of  water-cement  ratio  w/c  [l] : 

log10  Dq-  = 6.0  w/c  — 13.84  (1) 

Using  this  equation,  the  estimated  diffusivity  for  the  tank  is  3.6xl0-12  m2/s.  This, 
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along  with  a water-cement  ratio  of  0.40,  a 0.85  degree  of  hydration  (an  assumption  in 
the  original  4SIGHT  model),  and  a relation  for  the  porosity  0 [1], 


1 + 1.16a- 
1 + 3.2w/c 


(2) 


yields  an  estimate  for  the  capillary  porosity  of  0. 13,  which  is  consistent  with  the  measured 
values  of  porosity. 

The  4SIGHT  program  was  then  used  to  estimate  the  porosity.  The  program  started 
with  the  following  input  data  file  (See  Ref.  [1]  for  a description  of  each  parameter): 


DIFF  = 3.6B-12 
WC  = .40 

THICKNESS  = 0.100 
EXTERNAL  Na  = 0.00000 
EXTERNAL  K = 0.00000 
EXTERNAL  Cl  = 0.00000 
INTERNAL  K = 0.10000 
INTERNAL  Na  - 0.05000 
OUTPUT  K 
OUTPUT  Ca 
OUTPUT  Ca  OH 
REBAR  = .8000 
HEAD  =0.0 
DEPTH  = .2000 

JOINT  = 0.00000  AT  1.0  UNTIL  100. 
CRACK  = 0.00000  AT  2.0  DEPTH  0.25 
TIME  = 33000 


The  diffusivity  DIFF  is  in  units  of  m2/s;  the  concrete  thickness  THICKNESS  is  in  units  of 
meters;  the  INTERNAL  and  EXTERNAL  ion  concentrations  are  in  units  of  moles  per  liter. 

This  approach  assumes  that  there  were  no  cracks  or  joints  in  the  concrete  mortar. 
Also,  the  REBAR  and  DEPTH  variable  were  set  to  large  values  to  eliminate  the  effects  of 
corrosion  or  sulfate  attack. 

The  EXTERNAL  parameters  yield  a neutral  pH  external  environment: 


Initial  state  of  system: 


II 

ION 

" "EXTERNAL" 

"INTERNAL" 

II 

H 

" 0.00000 

0.00000 

II 

Ca 

" 0.00000 

0.00035 

II 

Na 

" 0.00000 

0.05000 

II 

K 

" 0.00000 

0.10000 

II 

OH 

" 0.00000 

0.15070 

II 

pH 

" 7.00000 

13.17810 

The  results  using  the  4SIGHT  estimate  for  the  diffusivity  follow: 
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W/C  =0.40 
D = 3.6  E -12 

Final  System  state: 


L (m) " 

"porosity' 

1 "pH" 

"fc" 

"Ca" 

"K" 

"CaOH 

0.0000 

0.1622 

7.000 

4525 

0.0000 

0.0000 

0.000 

0.0050 

0.1483 

12.404 

4747 

0.0123 

0.0003 

39.861 

0.0100 

0.1345 

12.411 

4990 

0.0119 

0.0009 

42.963 

0.0150 

0.1346 

12.421 

4989 

0.0114 

0.0017 

42.944 

0.0200 

0.1346 

12.431 

4988 

0.0109 

0.0024 

42.933 

Based  upon  the  values  for  the  porosity  and  the  calcium  hydroxide,  4SIGHT  predicts  that 
the  depth  of  leaching  is  between  5 mm  and  10  mm.  This  agrees  with  the  measured  values 
of  porosity.  However,  the  porosity  near  the  surface  is  only  a few  percent  greater  than  the 
bulk  value  of  13.5  %. 

A slight  modification  to  the  concrete  diffusivity  yields  a marked  change  in  the  4SIGHT 
estimate: 

w/c  = 0.40 
D = 5.0  E -12 


Final  System  state: 


"L(m) " 

"porosity"  "pH" 

"fc" 

"Ca" 

"K" 

"CaOH 

0.0000 

0.4123 

7.000 

2471 

0.0000 

0.0000 

0.000 

0.0050 

0.2896 

12.400 

3456 

0.0125 

0.0000 

9.814 

0.0100 

0.1670 

12.400 

4991 

0.0125 

0.0000 

31.971 

0.0150 

0.1671 

12.404 

4989 

0.0123 

0.0003 

31.953 

0.0200 

0.1672 

12.408 

4988 

0.0121 

0.0005 

31.942 

Increasing  the  concrete  diffusivity  from  3.6  to  5.0  (xl0~12  m2/s)  increased  the  estimated 
bulk  porosity  to  17  %;  since  the  correlation  between  porosity  and  diffusivity  has  con- 
siderable uncertainty,  it  is  difficult  to  assert  that  one  of  these  estimates  is  unfounded. 
The  depth  of  leaching  was  still  between  5 mm  and  10  mm.  However,  the  porosity  at 
the  surface  is  now  approximately  thrice  the  bulk  porosity,  as  was  seen  in  the  published 
measurements.  Although  a bulk  porosity  of  17  % is  greater  than  the  measured  value,  the 
size  limit  on  the  measured  capillary  pore  sizes  would  suggest  that  additional  capillary 
porosity  would  exist  in  pores  smaller  than  the  1 pm  limit  used  for  the  measurements.  At 
best,  the  values  of  diffusivity  appear  to  be  reasonable  upper  and  lower  bounds. 

2.1.2  Chloride  Diffusion 

The  Swedish  National  Testing  and  Research  Institute  initiated  a research  program  in 
which  concrete  specimens  were  placed  on  pontoons  in  sea  water  at  the  Traslovslage  Ma- 
rine Field  Station.  Concrete  specimens  from  more  than  50  mixture  designs  (0.25<wc<0.75) 
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Distance  (mm) 

Figure  2:  Chloride  concentration  (percent  by  mass  fraction  of  cement)  for  the  three 
concrete  mixtures  (solid  symbols).  The  concretes  are  made  with  ordinary  sulfate  resistant 
cements.  The  curves  are  estimates  by  4SIGHT.  Data  are  reproduced  from  Sandberg  et 
al.  [6]. 


were  arranged  in  three  tiers  that  define  their  exposure  zone:  submerged  zone,  splash  zone, 
and  atmospheric  zone.  A report  of  chloride  ingress  after  an  exposure  of  five  years  was  re- 
ported by  Sandberg,  Tang,  and  Andersen  [6].  Samples  were  extracted  from  the  exposure 
site  and  milled,  starting  from  the  exposure  surface,  in  1 mm  increments.  The  material 
at  each  depth  was  analyzed  for  total  chloride  content. 

The  report  lists  a number  of  concrete  mixtures  studied  in  the  program,  many  of  which 
contain  pozzolanic  additions.  Here,  only  the  results  for  the  ordinary  (sulfate  resistant) 
Portland  cement  (SRPC)  concrete  are  considered.  The  three  mixtures  are  denoted  SRPC 
0.3,  SRPC  0.4,  and  SRPC  0.5,  where  the  number  indicates  the  water-cement  ratio.  Also, 
since  4SIGHT  assumes  the  concrete  was  continuously  saturated,  only  the  data  for  the 
submerged  exposure  are  considered. 

In  the  report  of  Sandberg  et  al.  [6],  the  chloride  ion  diffusion  coefficient  was  estimated 
from  the  chloride  concentration  profile  within  the  concrete.  These  concentration  profiles 
are  shown  as  filled  symbols  in  Fig.  2.  It  is  common  that  chloride  concentrations  have  a 
maximum  within  the  concrete,  so  least  squares  regression  between  the  solution  to  Ficks 
law  and  the  recorded  values  was  applied  to  the  data  on  the  “falling”  side  of  the  curve. 
The  results  of  these  regression  analyses  are  shown  in  Table  1,  along  with  the  4SIGHT 
predictions  for  the  diffusion  coefficient  and  the  paste  porosity,  based  solely  upon  the 
water-cement  ratio.  The  curves  in  Fig.  2 are  the  4SIGHT  predictions  for  the  chloride 
concentration,  which  are  based  on  the  estimated  diffusion  coefficient,  not  the  reported 
value.  The  surface  chloride  concentration  for  the  4SIGHT  estimate  was  chosen  to  most 


Table  1:  Estimated  and  reported  diffusion  coefficients,  and  estimated  paste  porosities, 
for  the  three  concrete  mixtures  considered. 


4SIGHT 

Report 

Mixture 

DC\~ 

(10-12  m2/s) 

Porosity 

(paste) 

DC\~, 

(10  12  m2/s) 

SRPC  0.3 

0.9 

0.04 

1.5 

SRPC  0.4 

3.6 

0.14 

1.4 

SRPC  0.5 

14.5 

0.22 

2.3 

reasonably  approximate  the  data;  no  formal  regression  was  used.  Qualitatively,  only  the 
estimation  for  the  SRPC  0.3  mixture  is  accurate.  The  results  for  the  SRPC  0.4  mixture 
are  reasonable,  but  not  remarkable.  The  diffusion  coefficients  from  both  the  4SIGHT  and 
the  Sandberg  reports  are  shown  in  Table  1.  The  variation  in  the  diffusion  coefficients 
calculated  by  regression  analysis  [6]  was  less  than  expected.  This  is  striking  given  the 
range  in  estimated  paste  porosities. 

The  reason  for  the  similarity  among  diffusion  coefficients  calculated  by  Sandberg  et 
al.  [6]  can  be  explained  by  the  presence  of  magnesium  in  the  sea  water.  The  magnesium 
reacts  with  the  available  calcium  hydroxide  to  form  magnesium  hydroxide.  The  reaction 
products  consume  a greater  volume  than  the  reactants,  and  the  pores  become  filled  [7]. 
The  greater  the  quantity  of  calcium  hydroxide,  which  is  proportional  to  the  water-to- 
cement.  ratio,  the  greater  is  this  effect.  Since  a 0.3  water-cement  ratio  concrete  would 
have  relatively  little  available  calcium  hydroxide,  the  4SIGHT  estimation  for  this  mixture 
was  relatively  accurate.  This  effect  also  explains  why  the  reported  diffusion  coefficients 
for  the  other  concretes  were  much  lower  than  the  4SIGHT  estimates. 

2.1.3  HETEK  Model 

The  Hgjkvahtetsbeton-Entreprengrens  TEKnolgi  (HETEK)  {High  quality  concrete,  the 
Contractors  TEChnology}  program  was  initiated  by  the  Danish  Road  Directorate  to  es- 
tablish a system  for  estimating  concrete  service  life  for  exposure  to  marine  environments. 
The  entire  HETEK  program  is  subdivided  into  a number  of  tasks,  one  of  which  addresses 
chloride  penetration  into  concrete.  Although  HETEK  was  designed  primarily  for  con- 
crete in  marine  environments,  it  would  be  useful  to  simply  compare  the  HETEK  and  the 
4SIGHT  models. 

Since  the  HETEK  program  was  developed  for  a specific  purpose,  the  approach  to 
concrete  degradation  has  a focused  approach.  In  a marine  environment,  cements  with 
reduced  sulfate  content  are  specified  in  order  to  reduce  the  reaction  of  magnesium  from 
within  the  sea  water  with  sulfate  within  the  cement.  The  product  of  this  reaction,  mag- 
nesium sulfate,  leads  to  expansion.  An  additional  effect  due  to  the  low  sulfate  content  is 
that  the  magnesium  is  available  to  react  with  the  calcium  hydroxide  within  the  cement 
paste.  The  product  of  this  reaction,  magnesium  hydroxide,  has  a smaller  molar  volume 
than  calcium  hydroxide,  so  it  does  not  lead  to  expansion.  Rather,  the  magnesium  hy- 
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Table  2:  The  measured  diffusion  coefficients  (^jjpypiy)  and  the  estimated  diffusion 
coefficient  (F*4SIGHt)  f°r  different  values  of  water-cement  ratio  (vv/c). 


w/c 

dhetek 

D4SIGHT 

0.32 

3.9 

1.2 

0.40 

8.5 

3.0 

0.70 

25. 1 

229. 1 

droxide  seems  to  “clog"  the  pore  space,  effectively  reducing  the  diffusion  coefficient,  of 
the  concrete. 

Based  upon  experimental  data,  the  apparent  reduction  in  the  diffusion  coefficient  due 
to  the  magnesium  hydroxide  is  given  by  the  following  relationship  [2]: 


D(1)  = ( t0  \<> 

D{t0 ) vfex  + t) 


(3) 


The  quantity  tex  is  the  time  at  which  the  concrete  was  exposed  to  the  marine  envi- 
ronment. The  quantity  t()  is  the  time  at  which  the  concrete  diffusion  coefficient  has 
reached  a steady  value,  typically  6 months.  The  parameter  (3  varies  from  0.37  to  0.23  for 
water-cement  ratios  ranging  from  0.30  to  0.75,  respectively. 

The  HETEIv  model  is  based  upon  the  solution  of  Fields  equation  with  a time- 
dependent  diffusion  coefficient  D(t ): 


dc 

dt 


= D(t ) 


d2c 


dx2 


The  solution  to  this  equation  has  been  published  previously  [8]  for  diffusion  coefficients 
having  a time  dependence  as  given  in  Eqn.  3 above.  The  time  dependent  behavior 
for  a specific  mixture  is  based  upon  measurements  performed  on  similar  concretes  that 
have  had  similar  exposures.  The  result  is  a manual  [9]  that  consists,  in  large  part,  of 
nomographs  for  the  estimation  of  concrete  service  life,  defined  as  the  time  until  failure 
due  to  corrosion  of  the  steel  reinforcement. 

To  develop  the  analytical  equations  for  the  time-dependent  diffusion  coefficient,  the 
researchers  required  large  quantities  of  diffusion  data.  Only  a fraction  of  these  data  are 
shown  in  the  reports.  As  a basic  check,  a comparison  of  the  measured  diffusion  coefficients 
with  the  estimates  by  4SIGHT  based  upon  the  water-cement  ratio  is  given  in  Table  2.  The 
estimates  for  the  0.32  and  the  0.40  water-cement,  ratio  concretes  are  within  an  order  of 
magnitude  of  the  HETEIv  values.  This  is  encouraging  given  the  single  datum,  w/c,  used 
for  the  estimate.  The  4SIGHT  estimate  for  the  0.70  water-cement  ratio  concrete  differs  by 
an  order  of  magnitude.  Fortunately,  0.70  water-cement,  ratio  concretes  are  undesirable 
for  structures  designed  to  inhibit,  advective  and  diffusive  flow. 
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3 Transport 

Validating  the  transport  model  within  4SIGHT  is  paramount  to  validating  its  estimate 
of  service  life.  Every  degradation  mechanism  considered  by  4SIGHT  depends  on  the 
transport  of  ionic  species  through  the  concrete.  Therefore,  a proper  validation  requires 
that  the  transport  model  is  sufficiently  comprehensive  in  order  that  it  can  predict  the 
behavior  of  electrolytic  solutions. 

The  first  characteristic  of  electrolytic  solutions  to  be  considered  is  the  coulomb  in- 
teractions among  the  diffusing  ionic  species.  To  account  for  ion-ion  interactions  one 
must  determine  how  the  ions  would  respond  to  charge  imbalances.  The  solution  is  the 
electro-diffusion  equation,  which  has  been  used  extensively  in  solid  state  physics  [10],  the 
biological  sciences  [11,12],  and  only  recently  in  the  concrete  materials  research  [13-16]. 

The  second  characteristic  is  due  to  changes  in  the  solution  due  to  changes  in  species 
chemical  activity.  Since  the  pore  solution  of  cement  paste  typically  has  an  ionic  strength 
of  approximately  0.5  mol/kg,  one  must  also  account  for  the  nonideal  diffusive  behavior 
of  each  species.  This  nonideal  behavior  is  related  to  excluded  volume  and  the  second 
virial  coefficient.  This  effect  is  quantified  through  the  chemical  activity  of  the  species. 
The  chemical  activity  a of  a species  can  be  related  to  the  ionic  concentration  c through 
the  mathematical  construct  of  an  activity  coefficient  7: 

a = 7 c (4) 

The  value  of  the  activity  coefficient,  a dimensionless  number  between  zero  and  one, 
depends  on  the  concentration  and  units  used  for  the  concentration:  mol/L,  mol/kg, 
mol/m3,  etc.  Although  there  are  a number  of  approximation  schemes  for  calculating 
the  activity  coefficient,  the  Pitzer  equations  [17]  were  used  in  4SIGHT  since  relevant 
coefficients  are  readily  available  and  the  resulting  estimate  is  accurate  over  the  range  of 
ionic  strengths  typically  found  in  cement  paste  pore  solution. 

The  combination  of  the  electrostatic  interactions  and  the  activity  coefficient  calcula- 
tion yielded  a robust  transport  equation.  The  electro-diffusion  equation  can  be  used  to 
predict  the  single  diffusion  coefficient  for  a binary  salt,  which  are  reported  in  chemistry 
reference  handbooks,  over  a range  of  concentrations,  demonstrating  the  concentration 
dependence  on  the  apparent  diffusion  coefficient. 

3.1  Transport  Equation 

The  pore  solution  of  hardened  cement  paste  has  an  ionic  strength  of  approximately  0.5 
mol/kg.  At  this  ionic  strength,  the  pore  solution  cannot  be  modeled  using  equations  that 
characterize  idealized  transport,  which  generally  treat  each  ionic  species  independently  of 
one  another.  Further,  regardless  of  the  ionic  strength,  the  electrostatic  interactions  among 
the  various  species  are  quite  strong.  Although  the  self-diffusion  coefficient  for  sodium 
ions  is  smaller  than  that  of  chloride  ions,  the  electrostatic  interactions  enforce  electro- 
neutrality,  and  the  two  ions  diffuse  with  the  same  apparent  self-diffusion  coefficient,  giving 
a diffusion  coefficient  for  the  salt;  these  diffusion  coefficients  are  reported  in  chemistry 
reference  handbooks. 
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Further,  the  chemistry  handbooks  report  the  salt  diffusion  coefficient  as  a function 
of  concentration.  The  electrostatic  interactions,  by  themselves,  are  insufficient  to  char- 
acterize this  phenomenon.  One  approach  to  capture  this  behavior  in  a transport  model 
is  to  use  self-diffusion  coefficients  that  are  a function  of  concentration.  However,  the 
concentration  dependence  of  sodium  chloride  is  not  the  same  as  for  magnesium  chloride. 
Therefore,  some  experimenters  have  used  empirical  relations  to  predict  diffusion  coeffi- 
cients based  on  both  concentration  and  speciation  [18].  This  approach,  however,  can  only 
be  used  to  interpolate  results;  any  new  species,  or  concentrations  outside  the  range  from 
which  the  coefficients  were  determined,  would  invalidate  the  result. 

The  alternative  approach  is  to  characterize  the  ion-ion  interactions  at  a fundamental 
level,  thereby  “calculating”  the  effects  of  concentration  and  speciation.  The  electrostatic 
interactions  give  rise  to  an  electrical  held  that  ensures  electro-neutrality.  This  held  is 
unique  and  can  be  calculated  using  principles  of  electrostatics.  The  macroscopic  effect,  of 
the  electrostatic  interactions  is  an  electrical  potential  (the  diffusion  potential)  that  exists 
across  the  specimen,  and  is  measurable. 

The  concentration  dependence  arises  from  mechanical  ion-ion  interactions  that  only 
manifest  themselves  at  elevated  concentrations.  The  presence  of  a concentrated  ionic 
species  acts  to  change  the  hydrodynamic  properties  of  the  electrolyte,  thereby  changing 
the  apparent  viscosity  of  the  fluid.  This  behavior  is  related  to  the  species  activity,  which 
is  a function  of  its  chemical  potential. 

3.1.1  Electrical  Measurements 

Before  discussing  the  transport  equation,  it  will  be  useful  to  consider  the  relationship 
between  conductivity  and  diffusivity.  These  principles  are  fundamental  to  the  final  ex- 
pression of  the  transport  equation.  The  final  form  exploits  the  formation  factor,  which  is 
a material  parameter  that  can  be  determined  directly  from  conductivity  measurements. 

Conductivity  measurements  hold  great  promise  for  estimating  the  formation  factor 
of  saturated  porous  materials.  The  two  most  commonly  used  conduction  techniques  are 
the  electrical  migration  test  (driven  diffusion)  and  the  conductivity  measurement.  Mi- 
gration tests  are  used  to  determine  the  ionic  mobility,  which  can  be  related  to  diffusivity. 
Conduction  tests  are  used  to  estimate  the  formation  factor,  which  can  also  be  related  to 
diffusivity.  While  there  have  been  many  reports  on  the  use  of  the  electrical  migration 
tests  to  determine  the  diffusion  coefficient  of  cementitious  systems  [19],  there  have  been 
fewer  reports  on  the  use  of  the  formation  factor.  This  is  unfortunate,  since  the  formation 
factor  is  directly  related  to  the  diffusivity  of  the  material  pore  structure.  The  relationship 
between  formation  factor  and  diffusivity,  however,  has  subtleties  that  must  be  considered 
when  making  estimates  of  one  from  the  other. 

The  formation  factor  F has  its  origin  in  geological  research  on  saturated  porous  ma- 
terials. For  a nonconducting  porous  solid  saturated  with  a conducting  pore  solution,  the 
formation  factor  is  the  ratio  of  the  pore  solution  conductivity  ap  to  the  bulk  conductivity 
[20]: 

(5) 

This  quantity  characterizes  the  solid  microstructure  since  the  only  difference  between  the 
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conductivities  is  due  to  the  restricted  pathways  through  which  the  current  is  constrained 
in  the  bulk  conductivity  measurement. 

The  use  of  conduction  tests  to  estimate  the  diffusion  coefficient  is  due  to  the  relation- 
ship between  conductivity  and  diffusivity.  In  an  electrolytic  solution,  the  contribution  an 
ionic  species  makes  to  the  overall  conductivity  can  be  expressed  as  a function  of  its  con- 
ventional (electrochemical)  mobility  u,  amount-of-substance  concentration  c,  and  valence 
* [21]: 

a — z C-Tu  (6) 

The  quantity  T is  the  Faraday  constant.  The  mobility  can  then  be  related  to  diffusivity 
through  the  Einstein  relation  [21]: 


zTD  = RTu 


(7) 


The  quantity  R is  the  universal  gas  constant  and  T is  the  absolute  temperature.  Using 
these  equations,  one  can  estimate  the  diffusion  coefficient  from  either  measurements  of 
ion  mobility  (driven  diffusion)  or  from  the  conductivity  contribution  of  a particular  ion 
(conductivity). 

Alternatively,  if  a porous  material  is  saturated  with  a dilute  electrolytic  solution,  the 
ratio  of  the  pore  solution  conductivity  ap  to  the  bulk  conductivity  a b would  be  equal  to 
the  ratio  of  the  diffusion  coefficient  of  an  ion  in  the  pore  solution  Dp  to  the  bulk  diffusion 
coefficient  Db  of  that  ion: 


F 


Op 


Dp 

Db 


(8) 


Note  that  this  relationship  holds  even  if  the  Einstein  relation  is  incorrect  at  large  con- 
centrations; the  multiplicative  error  occurring  in  both  Dp  and  Db  would  cancel.  Since 
the  self-diffusion  coefficient  of  ions  can  be  found  in  books,  this  is  a provocative  approach 
to  determining  the  bulk  diffusion  coefficient  Db  [22]. 

Unfortunately,  the  pore  solution  of  typical  cementitious  systems  can  have  a large  ionic 
strength  (0.5  mol-kg-1  to  1.0  mol-kg-1),  and  so  care  must  be  exercised  when  estimating 
the  bnlk  diffusion  coefficient  from  the  formation  factor.  At  these  large  ionic  strengths, 
the  pore  solution  diffusion  coefficient  Dp  of  an  ion  is  different  from  the  self-diffusion 
coefficient,  which  is  the  value  reported  in  tables  [23].  Equation  8 can  still  be  used,  if  one 
can  correctly  determine  the  diffusion  coefficient  Dp  in  the  pore  solution.  Unfortunately, 
given  the  constraints  of  most  experimental  diffusion  apparatus,  and  that  experimenters 
typically  measure  time  dependent  concentrations,  an  exact  value  for  Dp  can  be  difficult 
to  define. 

An  alternative  approach  is  to  determine  the  formation  factor  from  the  diffusion  data 
by  separating  the  microstructural  dependence  from  the  concentration  dependence  in  the 
diffusion  coefficient.  It  will  be  shown  that  the  formation  factor  can  be  alternatively  ex- 
pressed as  a function  of  the  microstructural  diffusion  coefficient  of  the  porous  material. 
The  microstructural  diffusion  coefficient  is  independent  of  both  changes  in  the  pore  so- 
lution chemistry  and  the  diffusion  potential.  The  estimated  formation  factor  could  then 
be  used  to  estimate  the  diffusive  transport  of  any  ion  within  the  studied  system  since  it 
uniquely  characterizes  the  solid  pore  structure:  the  chemical  effects  would  be  calculated 
independently. 
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3.1.2  Diffusive  Ionic  Flux 


The  proper  application  of  the  formation  factor  to  porous  systems  containing  concentrated 
electrolytes  can  be  best  studied  by  starting  with  the  appropriate  transport  equation. 
For  the  case  of  diffusing  charged  species  in  concentrated  electrolytes,  the  transport  is 
governed  bv  coupled  electrical  and  diffusive  transport,  and  ion-ion  interactions  must  be 
also  be  considered. 

The  diffusion  of  ionic  species  in  an  electrolyte  is  governed  by  the  electro-diffusion 
equation.  For  the  z-th  ionic  species,  the  electro-diffusion  equation  relates  the  bulk  flux 
j to  the  concentration  c,  the  diffusion  potential  ipD,  the  bulk  microstructural  diffusion 
coefficient  D^,  and  the  bulk  conventional  mobility  u [10]: 

j - -Dfj.  ^1  + Vc  ~ zTucViPd  (9) 

The  subscript  z has  been  omitted  here  in  order  to  improve  visual  clarity  of  the  equation, 
and  will  be  added  later.  The  quantity  7 is  the  activity  coefficient  for  the  species.  Al- 
though this  equation  neglects  adsorption  effects,  it  is  otherwise  a complete  description 
of  the  nonreacting  diffusive  transport  of  charged  species  in  concentrated  electrolytes. 
Equation  9 still  bears  a resemblance  to  Fick’s  law  for  apparent  diffusion  Da  [24]: 

j = -DaVc  (10) 

The  quantity  D^{  1 TcHny/cMnc]  is  an  agglomerated  diffusion  coefficient  that  includes 
the  effects  of  changes  in  both  the  microstructure  and  the  concentration. 

It  should  be  noted  that  the  agglomerated  diffusion  coefficient  is  not  the  apparent 
bulk  diffusion  coefficient  D Strictly  speaking,  the  apparent  diffusion  coefficient  also 
includes  the  effects  of  the  diffusion  potential  i/jp.  In  those  cases  where  the  self-diffusion 
coefficients  of  all  the  diffusing  species  are  nearly  identical,  the  diffusion  potential  will  be 
nearly  zero,  and  the  apparent  diffusion  coefficient  will  be  nearly  equal  to  the  agglomerated 
diffusion  coefficient.  However,  for  cementitious  systems,  there  are  many  species  present 
with  varying  self-diffusion  coefficients. 

Ideally,  one  wonld  like  to  distinguish  between  the  effects  due  to  microstructural 
changes  and  the  effects  due  to  changes  in  the  pore  solution  electro-chemistry.  By  ob- 
servation, one  can  see  that  the  microstructural  diffusion  coefficient  D M characterizes  the 
solid  microstructure  and  that  the  quantity  in  parenthesis  characterizes  the  pore  solution 
chemistry.  The  microstructural  diffusion  coefficient  is  itself  independent  of  the  pore  solu- 
tion chemistry.  A demonstration  of  how  to  extract  the  value  of  microstructural  diffusion 
coefficient  from  a diffusion  experiment  is  accomplished  through  a detailed  discussion  of 
the  formation  factor. 

3.1.3  Formation  Factor 

Consider  an  electrical  measurement  of  the  formation  factor  F performed  on  a porous 
specimen  saturated  with  a concentrated  electrolyte.  The  conductivity  of  the  pore  solution 
Op  is  a function  of  the  ionic  strength  /.  The  conductivity  is  proportional  to  the  ionic 
strength  such  that  as  the  ionic  strength  approaches  a value  of  zero,  the  conductivity  also 
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approaches  a value  of  zero.  Similarly,  the  hulk  conductivity  has  the  same  dependence 
on  the  pore  solution  conductivity.  If  the  pore  solution  conductivity  doubles,  the  bulk 
conductivity  will  also  double,  ignoring  surface  conduction  contributions. 

The  diffusion  coefficient  of  an  ion  in  a concentrated  electrolyte  has  a different  type 
of  dependence  on  the  ionic  strength.  The  self-diffusion  coefficients  reported  in  tables  are 
limiting  values  under  the  condition  that  the  limit  that  both  the  concentration  and  the 
diffusion  potential  go  to  zero.  In  the  practical  dilute  limit,  the  cations  and  the  anions  of 
a 1:1  valence  electrolyte  must  diffuse  together  due  to  charge  neutrality.  This  leads  to  the 
creation  of  the  diffusion  potential.  Therefore,  the  pore  solution  diffusion  coefficient  Dp 
must  be  an  function  of  the  ionic  strength,  act  ivity  coefficients,  and  the  diffusion  potential: 

Dp  = Ds  [1  + g{F  7,  -0d)]  (11) 

The  function  g is  an  arbitrarily  complex  function  of  the  ionic  strength  /,  the  activity 
coefficient  7,  and  the  diffusion  potential  Uz>-  It  characterizes  the  interaction  between 
diffusion  ionic  species  and  the  bulk  electrolyte.  The  quantity  Ds  is  the  self-diffusion 
coefficient  reported  in  tables  [23].  A similar  relationship  characterizes  the  bulk  diffusion 
coefficient  Db : 

Db  = Dp  [1  + g(I,  7,  i/)d)\  (12) 

The  function  g is  the  same  since  the  diffusion  ionic  species  will  see  the  same  electrolyte. 
The  microstructural  diffusivity  Dp  is  the  value  corresponding  to  the  dilute  limit  for 
the  bulk  diffusivity.  It  is  independent  of  both  the  activity  coefficient  and  the  diffusion 
potential,  and  is  the  same  quantity  that  appears  in  Eqn.  9. 

Using  these  relations,  the  formation  factor  can  be  written  as  a function  of  quantities 
that  are  themselves  functions  of  the  ionic  strength  /,  activity  coefficient  7,  and  the 
diffusion  potential  tbo- 

= <7P(I)  = Ds  [1  + g(F 7,  M]  = Ds  ,13v 

ab(I)  Dp  [1  + g(1 , 7,tM]  Dn 

However,  the  resulting  formation  factor  F is  independent  of  both  the  ionic  strength  and 
the  ionic  species;  a fact  that  is  exploited  by  numerical  calculations  of  formation  factor 
on  model  microstructures  [25]. 

The  advantage  of  using  the  ratio  Ds/Dp  is  not  immediately  obvious.  In  Eqn.  8,  Db 
is  observable,  but  Dp  is  difficult  to  determine  precisely  in  typical  diffusion  experiments 
containing  a variety  of  species,  each  with  varying  concentration.  In  Eqn.  13,  Ds  can 
be  found  in  tables,  but  the  quantity  Dp  cannot  be  observed  in  electrolytes  because  the 
effects  of  the  diffusion  potential  cannot  be  eliminated.  Serendipitously,  the  usefulness  of 
Dp  in  service  life  prediction  through  its  use  in  the  electro-diffusion  equation  is  also  the 
key  to  determining  this  quantity  from  experiment. 

The  format  ion  factor  is  extremely  useful  in  models  of  ionic  transport  . Accurate  service 
life  predictions  require  models  that  can  account  for  independent  changes  in  either  the  solid 
microstructure  or  the  pore  solution  chemistry.  A suitable  computer  program  could  solve 
the  electro-diffusion  transport  equation  by  using  estimates  of  the  ion  activity  coefficient  7 
from  published  empirical  relations  and  using  the  local  electro-neutrality  condition  to  solve 
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for  the  diffusion  potential.  The  microstructural  diffusion  coefficient  in  the  electro- 
diffusion  equation  would  be  calculated  from  the  formation  factor  F and  the  dilute  limit 
self  diffusion  coefficient  Ds: 

Ds 

D,  = t (14) 

Note  that  this  relation  holds  for  all  the  species,  yet  there  is  only  a single  value  for  the 
formation  factor  F. 

Alternatively,  one  can  determine  the  formation  factor  from  experimental  diffusion 
data  by  using  the  same  computer  program.  Using  a computer  program  that  implements 
Eqns.  9 and  14,  one  adjusts  the  formation  factor  until  the  computed  output  matches  the 
experimentally  observed  quantities.  The  experimental  program  described  subsequently 
used  this  approach  to  determine  the  formation  factor  from  divided  cell  diffusion  mea- 
surements. 

3.1.4  Conservation  Equation 


p, 


Specimen 


v, 


v 


P 


Figure  3:  Schematic  of  the  bulk  flux  and  the  pore  flux,  characterized  by  the  bulk  velocity 
vb  and  pore  velocity  vp,  due  to  a pressure  hydrostatic  pressure  drop  of  P\  — P2  across 
the  specimen. 

The  flux  in  Eqn.  9 was  developed  for  transport  within  an  electrolyte.  Transport 
through  a porous  material  must  account  for  the  solid  framework  that  comprises  the 
material.  The  schematic  of  the  Darcy  experiment  shown  in  Fig.  3 demonstrates  the 
difference  between  the  pore  flux  (through  electrolyte)  and  bulk  flux  (through  porous 
material). 

The  experiment  shown  in  Fig.  3 includes  a porous  specimen  saturated  with  a brine. 
The  volumes  above  and  below  the  specimen  are  filled  with  the  same  brine.  Due  to 
a hydrostatic  pressure  Pi  applied  above  the  sample,  there  is  a mass  flow  of  the  brine 
through  the  porous  specimen.  The  top  free  surface  of  fluid  “falls”  at  a velocity  of  vb, 


representing  the  bulk  velocity  of  the  fluid  through  the  sample.  If  the  fluid  above  the 
sample  were  instantaneously  filled  with  a colored  dye,  the  “front”  of  the  dye  would 
advance  through  the  specimen  at  the  pore  velocity  vp.  The  pore  flux  jp  of  ionic  species 
within  the  brine  would  be  proportional  to  the  species  concentration  c:  jp  = cvp.  Similarly 
for  the  bulk  flux:  }b  = cvb. 

Employing  conservation  of  mass,  and  assuming  that  water  is  incompressible  at  these 
pressures,  the  pore  velocity  vp  will  greater  than  the  bulk  velocity  vb.  The  ratio  of  the 
velocities  is  the  connected  porosity  6:  v6  = </>vp.  The  porosity  is  also  the  ratio  of  the  pore 
area  and  the  bulk  area:  Ap  = 4>Ab.  Likewise,  by  conservation  of  ionic  species,  the  rate  at 
which  species  are  transported  through  the  system  is  same  above,  inside,  and  below  the 
specimen: 


j b Ab  jp  Ap 

j/>  — 4*3  p 


(15) 

(16) 


This  is  an  important  relationship  when  considering  diffusive  and  advective  transport 
through  porous  media.  With  a relationship  between  the  pore  flux  jp  and  the  bulk  flux 
jft,  one  can  now  develop  the  appropriate  flux  relationships  for  each  of  the  three  transport 
mechanisms  (diffusion,  permeation,  and  conduction)  considered  in  4SIGHT. 

As  mentioned  previously,  the  hydrodynamic  flux  will  occur  due  to  Darcy  flow  under 
a hydraulic  pressure  gradient  VP  across  the  specimen.  The  bulk  velocity  vb  is  equivalent 
to  the  Darcy  velocity: 

Vt  = — VP  (17) 

A 

This  result,  and  the  similar  one  for  the  pore  velocity,  can  be  substituted  into  the  appro- 
priate flux  equation: 


3b 

^kb  „ 

= cvft  = V P 

A 

(18) 

jp 

cA/p 

= cvp  = - VP 

A 

(19) 

kb 

II 

?r- 

72 

(20) 

In  nearly  every  application  of  Darcy’s  law  to  concrete,  researchers  report  the  bulk  per- 
meability jfcft,  so  the  bulk  flux  jft  is  most  convenient. 

Although  it  has  been  argued  that  the  electro-diffusion  equation  is  a more  complete 
description  of  diffusive  transport,  the  simplified  Fick’s  law  is  used  here  as  a demonstration 
of  the  distinction  between  pore  flux  and  bulk  flux: 


3b  — 

— DbX7c 

(21) 

jp  ~ 

— DpVc 

(22) 

Db  = 

<t>Dp 

(23) 

Reported  results  from  diffusion  measurements  vary  between  reporting  either  the  pore 
flux  jp  or  the  bulk  flux  j6.  Unfortunately,  this  distinction  is  not  always  made,  and  so  the 
reader  must  exercise  care  when  interpreting  the  results  from  diffusion  experiments. 
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Figure  4:  Schematic  for  ionic  flux  conservation  through  a volume  of  unit  area  and  differ- 
ential width  dx. 


The  conductive  flux  is  that  which  occurs  due  to  the  application  of  an  electrical  field. 
Strictly  speaking,  the  flux  is  proportional  to  the  conventional  ionic  mobility  u.  The  pore 
flux  and  bulk  flux  depend  on  the  corresponding  pore  and  bulk  mobilities: 


it,  = -zcTubVi/j 

(24) 

jp  = — zctFupX/ ^ 

(25) 

Uh  = </>?ip 

(26) 

In  general,  the  mobilities  are  approximated  using  the  Einstein  relation: 


u — 


D 

RT 


The  diffusion  coefficient  D is  the  self-diffusion  coefficient  reported  in  tables  [23].  This 
relationship,  however,  is  only  valid  in  the  dilute  limit. 

With  the  pore  and  bulk  flux  equations  established  for  both  diffusion  and  conduction 
transport,  a time-dependent  equation  is  needed  in  order  to  calculate  changes  in  concen- 
tration. This  equation  can  be  derived  from  the  conservation  of  species.  Figure  4 shows  a 
volume  element  with  unit  area  and  differential  width  dx.  The  ^-component  of  the  species 
flux  entering  the  volume  is  jx(x,t ),  and  the  flux  exiting  is  jx(x  + dxj).  The  rate,  per 
unit  area,  at  which  species  are  filling  the  volume  is  jx{xj)  — jT(x  + dxJ).  The  change 
in  the  concentration  is  this  difference  divided  by  the  width  dx: 


dc  _ jx{x,t)  - j.T{x  + dx,t)  _ -djx 
dt  dx  dx 


dc 

dt 


-Vj 


(28) 


Now  this  equation  must  be  reformulated  to  account  for  a porous  material. 

The  equation  for  a porous  material  can  be  obtained  from  the  definition  of  concen- 
tration. The  experimentally  useful  quantity  is  the  intrinsic,  or  pore  concentration  c,  in 
units  of  moles  of  species  per  unit  pore  volume  Ip.  Let  the  bulk  concentration  c represent 
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the  number  moles  per  unit  bulk  volume  1 b (volume  of  pore  solution  and  solid  structure 
combined).  Apply  Gauss’s  divergence  theorem  to  Eqn.  28  over  the  bulk  volume  \b: 


/^dvi  = -/v'j‘du  (29) 

The  bulk  concentration  c can  be  related  to  the  intrinsic  concentration  through  the  poros- 
ity; c = (j)C.  The  resulting  conservation  equation  relates  the  bulk  flux  j/,  to  the  intrinsic 
concentration  c : 


d(f)c 

dt 

-V  j6 

(30) 

II 

-V  • 0jp 

(31) 

These  relationships  have  been  published  elsewhere  [26-29],  primarily  in  the  field  of  geo- 
physical research. 


3.1.5  Complete  Transport  Equation 

The  complete  transport  equation  has  separate  terms  for  the  diffusion,  conduction,  and 
permeation  components  of  the  flux.  Substituting  for  the  bulk  microstructural  diffusion 
D/x,  the  self-diffusion  coefficient  the  formation  factor  F,  and  the  definition  of  the 
Darcy  velocity  yields  the  following  time-dependent  transport  equation. 


d<f>Cj 

dt 


■V  ■ jfcii 


Vc‘  + IFF  IF  V ' c'V</’b  + V • ~VP  (32; 
HI  t n 


This  transport  equation  is  significant  in  that  only  three  numbers  (A,  F,  kb)  are  required 
in  order  to  characterize  the  transport  of  an  arbitrary  number  of  ionic  species. 
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4 Transport  Validation 

Validation  of  the  transport  equation  within  4SIGHT  is  a two-stage  process.  First,  4SIGHT 
will  be  used  to  calculate  the  properties  of  bulk  fluids.  These  properties  include  density, 
activity  coefficient,  and  salt,  diffusion  coefficient.  These  results  can  then  be  compared 
to  values  reported  in  chemistry  handbooks.  Second,  4SIGHT  will  be  used  to  calculate 
the  response  of  a controlled  diffusion  experiment  performed  on  a well-defined  porous 
material.  The  porous  material  is  nonreactive  and  has  relatively  little  surface  charge. 

4.1  Reference  Data 

Comparisons  between  4SIGHT  and  reference  data  are  needed  in  order  to  ensure  that  the 
fundamental  building  blocks  of  4SIGHT  are  constructed  in  an  orderly  manner.  Practical 
field  data  are  nearly  always  reported  in  units  of  molarity  (mol/L),  or  similar  units  of 
quantity  per  volume.  The  calculation  of  the  activity  coefficient,  however,  is  based  on 
concentration  expressed  in  units  of  molality  (mol/kg).  Converting  between  molarity 
and  molality,  requires  an  estimate  of  the  solution  density.  From  a coherent  means  of 
calculating  density  and  activity  coefficient,  one  can  then  use  the  electro-diffusion  equation 
to  calculate  the  diffusion  coefficient  of  various  binary  salts  and  compare  the  results  to 
values  reported  in  chemistry  handbooks. 

4.1.1  Electrolyte  Density 


Figure  5:  Density  of  1-1  electrolytes  (a)  and  1-2,  2-1,  and  2-2  electrolytes  (b)  as  a function 
of  molarity.  Data  from  the  chemistry  handbook  appear  as  symbols.  The  predictions  by 
4SIGHT  appear  as  solid  curves. 

A review  of  approximation  schemes  for  various  electrolyte  properties  revealed  no 
established  means  of  estimating  the  electrolyte  density.  The  scheme  implemented  by 
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4SIGHT  was  to  assume  that  each  ionic  species  contributes  independently  to  the  overall 
density.  The  contribution  is  a function  of  the  density  and  the  mole  fraction  of  that  species 
within  the  electrolyte. 

In  Figs.  5a  and  5b  are  the  results  for  nominally  pH  neutral  salts.  The  density  is 
shown  as  a function  of  electrolyte  concentration.  The  estimates  of  4SIGHT  are  denoted 
bv  solid  curves.  The  data  from  the  CRC  Handbook  [30]  are  denoted  by  Filed  symbols. 
The  results  for  the  1-1  electrolytes  are  shown  in  Fig.  5(a),  and  those  for  the  1-2,  2-1,  and 
2-2  electrolytes  are  shown  in  Fig.  5(b).  In  both  cases,  the  approximations  are  relatively 
accurate  over  the  full  range  of  concentrations.  In  general,  the  accuracy  for  the  1:1  salts 
are  better  than  for  the  1-2  and  2-1  salts.  For  both  cases,  the  accuracy  is  well  within  one 
percent  over  the  concentration  range  shown. 


Figure  6:  Density  of  strong  acids  and  bases.  Data  from  the  chemistry  handbook  appear 
as  symbols.  The  predictions  by  4SIGHT  appear  as  solid  curves. 

The  density  estimates  for  strong  acids  and  bases  are  shown  in  Fig.  6 as  a function 
of  electrolyte  concentration.  The  estimates  by  4SIGHT  are  denoted  by  solid  curves.  The 
handbook  values  for  acids  are  denoted  by  Filed  symbols,  and  values  for  bases  are  denoted 
by  open  symbols.  Generally,  the  results  for  acids  and  bases  are  not  as  accurate  as  for 
the  pH  neutral  salts.  Nevertheless,  the  accuracy  is  still  less  than  a percent  over  the 
concentration  range  shown. 

4.1.2  Activity  Coefficient 

Accurate  estimates  of  activity  coefficients  are  necessary  for  both  calculations  of  transport 
Fux  at  elevated  concentrations  and  for  calculations  of  equilibrium  species  concentrations 
in  reacting  systems.  In  4SIGHT.  the  activity  coefficient  is  calculated  using  the  equations 
of  Pitzer  [17].  To  reduce  complexity  and  gain  computational  speed,  a “truncated”  im- 
plementation was  used.  All  interactions  composed  of  two  and  three  charged  species  were 
implemented;  all  other  interactions  were  been  omitted. 
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Figure  7:  Activity  coefficient  of  chloride  ions  7 ci-  as  a function  of  molality  m for  various 
electrolytes.  The  estimates  by  4SIGHT  are  shown  as  solid  curves.  Results  from  the 
PHRQPITZ  computer  program  [31]  are  shown  as  filled  symbols.  Data  from  the  CRC 
Handbook  of  Chemistry  and  Physics  [30]  are  shown  as  open  symbols. 


As  a measure  of  performance,  the  results  from  4SIGHT  are  compared  to  both  CRC 
Handbook  data  and  the  results  of  the  PHRQPITZ  [31]  computer  program,  which  is  a full 
implementation.  Shown  in  Fig.  7 are  the  results  for  the  activity  coefficient  of  chloride  ions 
7 ci-  as  a function  of  electrolyte  concentration.  The  estimates  by  4SIGHT  are  denoted 
by  solid  curves.  The  values  from  PHRQPITZ  are  denoted  by  filled  symbols,  and  the 
values  from  the  CRC  Handbook  are  denoted  by  open  symbols.  Generally,  the  agreement 
among  all  three  is  very  good  over  the  entire  range  of  electrolyte  concentration  shown, 
suggesting  that  the  4SIGHT  implementation  is  sufficient.  This  example  of  chloride  ion 
activity  is  characteristic  of  the  accuracy  with  which  4SIGHT  can  estimate  the  activity  of 
all  the  ionic  species  incorporated  into  the  program. 

4.1.3  Electrolyte  Diffusion  Coefficient 

The  electrolyte  diffusion  coefficient  was  calculated  by  4SIGHT  and  compared  to  values 
published  in  the  CRC  Handbook.  The  “sample”  had  a porosity  of  one,  and  a formation 
factor  of  one.  The  change  in  concentration  across  the  specimen  was  only  a few  percent. 
The  flux  of  salt  passing  through  the  solution  was  then  used  to  calculate  an  apparent 
diffusion  coefficient,  based  on  Fick's  law  (Eqn.  10).  This  apparent  diffusion  coefficient 
was  then  compared  with  the  handbook  value. 

Data  for  various  electrolytic  solutions  are  shown  in  Table  3 over  varying  concentra- 
tions. Generally,  4SIGHT  is  accurate  to  within  a few  percent  up  to  0.1  mol/L.  The  error 
increases  with  concentration  and  is  greater  for  larger  valence  species.  The  results  are  still 
reasonable,  considering  that  a species  with  a 1.0  mol/L  concentration  within  concrete 
pore  solution  is  relatively  rare. 


21 


Table  3:  Comparison  of  diffusion  coefficients  D from  the  computer  program  (CP)  and 
handbook  (HB)  (Ref.  [30])  values  for  various  salts  in  aqueous  solution. 


Salt 

cone. 

mol-L-1 

Dhb 

10~9  m2-s_1 

Dcp 

10-9  m2-s-1 

Dcp-Oh  B 

Dhb 

KC1 

0.01 

1.917 

1.902 

-0.0078 

0.10 

1.844 

1.807 

-0.0201 

1.00 

1.892 

1.801 

-0.0481 

NaCl 

0.01 

1.545 

1.539 

-0.0039 

0.10 

1.483 

1.476 

-0.0047 

1.00 

1.484 

1.571 

+0.0586 

LiCl 

0.01 

1.312 

1.307 

-0.0038 

0.10 

1.269 

1.277 

+0.0063 

1.00 

1.302 

1.555 

+0.1943 

LiBr 

0.10 

1.279 

1.286 

+0.0055 

1.00 

1.404 

1.647 

+0.1731 

LiNCE 

0.01 

1.276 

1.277 

+0.0008 

0.10 

1.240 

1.242 

+0.0016 

1.00 

1.293 

1.467 

+0.1346 

kno3 

0.01 

1.846 

1.831 

-0.0081 

KI 

0.10 

1.865 

1.829 

-0.0193 

1.00 

2.065 

1.911 

-0.0746 

Na2S04 

0.01 

1.123 

1.059 

-0.0570 

MgS04 

0.01 

0.710 

0.568 

-0.2000 

The  result  for  MgS04  is  rather  anomalous.  The  value  calculated  by  4SIGHT  is  consis- 
tent with  the  approximation  based  on  the  Nernst  relation: 


D = (z+  + z_) 


D+D . 

z^D-^  -(-  Z—D _ 


d In  7±  \ 
8 Inc  ) 


(33) 


The  values  calculated  by  4SIGHT  are  consistent  with  this  equation,  as  are  many  of  the 
handbook  values.  The  discrepancy  in  the  MgS04  result  may  be  due  to  either  experi- 
mental errors  in  the  handbook  value  or  the  use  of  the  Einstein  relation  in  Eqn.  7 as  an 
approximation  for  the  ionic  mobility  u.  A more  rigorous  estimate  would  require  a con- 
certed effort  to  develop  an  estimate  of  the  ionic  mobility  as  a function  of  concentration. 


4,2  Laboratory  Data 

The  laboratory  experiments  were  performed  in  two  stages.  The  first  stage  validates  the 
use  of  the  formation  factor  in  the  transport  equation.  Independent  electrical  measure- 
ments of  formation  factor  are  compared  with  values  estimated  from  diffusion  data.  The 
second  stage  consists  of  a series  of  diffusion  measurements  that  attempts  to  demonstrate 
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that  the  transport  equation  is  valid  for  a range  of  ionic  species.  Further,  the  varied 
behavior  of  these  electrolytes  was  used  to  demonstrate  that  the  apparent  bulk  diffusion 
coefficient  depends  strongly  on  the  species  present,  and  comparatively  weakly  on  the 
ionic  strength,  over  the  range  of  concentrations  that  are  of  interest. 

Results  from  the  two  stages  of  experiments  have  been  presented  [32,33]  and  submitted 
for  publication  [34,35].  The  experimental  details  and  the  results  are  repeated  here  for 
completeness. 

4.2.1  Ceramic  Frit 

The  diffusion  experiments  were  performed  on  a commercial  porous  frit  composed  of  a 
sintered  alumina  ceramic  (99.8  % alumina)  measuring  (50.76T0.04)  mm  in  diameter  and 
(6.48T0.03)  mm  thick;  the  uncertainty  is  an  estimated  standard  deviation  calculated  from 
7 frits.  Mercury  intrusion  porosirnetry  analysis  revealed  a porosity  of  approximately  0.26, 
and  a pore  size  distribution  concentrated  between  350  nm  and  450  nm.  This  is  consistent 
with  the  manufacturer’s  published  specification  of  a pore  size  less  than  500  nm. 

The  alumina  frit  was  regarded  as  a nearly  ideal  porous  medium  for  studying  diffusive 
transport.  The  alumina  is  chemically  stable  at  moderate  values  of  pH.  in  addition, 
conductivity  measurements  suggest  that  the  frit  has  a relatively  low  bound  surface  charge. 
As  defined  in  Eqn.  5,  the  formation  factor  F is  calculated  from  the  ratio  of  the  pore  fluid 
electrical  conductivity  ap  to  the  the  bulk  electrical  conductivity  ab  with  that  pore  fluid. 
The  calculation  can  be  repeated  for  different  pore  solutions  with  varying  values  of  ap. 


0 10  20  30  40  50 
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Figure  8:  The  pore  solution  conductivity  op  and  the  corresponding  bulk  conductivity  ob 
for  the  resin  system  of  McDuff  and  Ellis  [36]  and  the  frit  system  used  here. 

Such  an  experiment  was  performed  in  1979  by  McDuff  and  Ellis  [36]  on  a model  porous 
material  composed  of  a packed  column  of  cation  exchange  resin.  Their  data  for  large 
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pore  solution  conductivities  are  denoted  in  Fig.  8 by  Filed  circles.  Linear  ordinary  least 
squares  regression  applied  to  the  data  reveal  a positive  intercept  of  (0.G8T0.17)  S/m;  the 
uncertainty  is  an  estimated  standard  deviation.  Similarly,  one  of  the  frits  was  saturated 
with  standard  potassium  chloride  conductivity  solutions  [37],  and  the  data  are  also  shown 
in  Fig.  8 as  Filed  triangles.  For  the  frit,  the  intercept  was  only  (0.00258T0. 00025)  S/m. 
Because  the  intercept  is  proportional  to  the  surface  charge  of  the  material  [38].  one 
expects  that  the  surface  charge  of  the  frit  is  approximately  250  times  smaller  than  that 
of  the  resin.  Unfortunately,  comparisons  with  cementitious  systems  are  not  possible  due 
to  lack  of  similar  data.  Regardless,  the  data  presented  herein  for  the  frit  system  will 
show  that  an  analysis  that  neglects  binding  still  captures  all  the  salient  features  of  the 
effect  that  speciation  has  on  the  diffusive  transport  of  ions  through  these  ceramic  frits. 

4.2.2  Experimental  Setup 

The  divided  cell  arrangement  was  used  as  the  experimental  apparatus  because  of  its 
procedural  simplicity.  Since  the  ceramic  specimens  have  a relatively  large  diffusion  co- 
efficient, compared  to  to  hydrated  cement  paste,  usable  data  could  be  collected  in  a 
reasonable  period  of  time.  As  mentioned  previously,  the  ceramic  material  also  has  a 
relatively  small  surface  charge  and  surface  adsorption  effects  could  be  neglected.  The 
result  is  a material  that  should  exhibit  nearly  ideal  diffusive  transport  properties. 

The  target  ion  for  the  experiment  was  iodide.  T he  iodide  concentration  was  deter- 
mined using  a commercial  solid-state  ion  selective  combination  electrode.  Iodide  was 
chosen  as  the  target  ion  because  of  experimental  advantages  iodide  electrodes  have  over 
chloride  electrodes.  For  each  measurement,  the  electrode  was  standardized  using  refer- 
ence iodide  solutions. 

Each  ceramic  frit  was  Frst  mounted,  using  epoxy,  into  an  acrylic  annulus.  After  the 
epoxy  had  cured,  the  specimen  was  saturated  with  a potassium  iodide  solution  under 
vacuum.  The  saturated  frit  was  then  clamped  between  two  glass  vessels,  each  having  a 
capacity  of  approximately  250  mL,  using  rubber  o-rings  for  a water  tight  seal.  One  vessel 
was  Filed  wit  h KI  and  the  other  with  a test  electrolyte.  Both  vessels  were  sealed  to  reduce 
evaporation,  and  then  placed  inside  a walk-in  environmental  chamber  maintained  at 
25  °C.  All  sampling  and  concentration  measurements  were  performed  within  the  chamber 
in  order  to  ensure  stability  in  both  the  diffusion  and  the  measurements;  there  is  a nominal 
2 %/°C  variation  in  the  self-diffusion  coefficient  [23]. 

4.2.3  Analysis 

Historically,  ionic  transport  in  cementitious  systems  has  been  characterized  by  Fick’s  law 
of  diffusive  transport: 

j,  = -Dye,  (34) 

The  quantity  represents  an  apparent  diffusion  coefficient  for  the  2-th  species.  The 
analysis  of  the  experimental  diffusion  data  is  based  on  the  “constant  gradient”  assumption 
that  underlies  Fick’s  law.  Based  on  this  assumption,  the  concentration  gradient  across 
the  specimen  is  constant  and  equal  to  the  difference  between  cell  concentrations,  divided 
by  the  specimen  thickness  L.  Under  these  conditions,  Fick’s  law  applies.  For  vessel 
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Figure  9:  Schematic  of  the  constant  gradient  hypothesis  across  a sample  with  thickness 
L.  Vessels  1 and  2 contain  an  aqueous  electrolyte  with  target  species  concentrations  C] 
and  c2,  and  have  volumes  V\  and  u2,  respectively. 


volumes  V\  and  u2  and  specimen  area  A , the  difference  in  concentration  A between  the 
two  vessels  can  be  shown  to  decay  exponentially  [35,39]: 


A 

An 


- exp 


-DM 

L 


(V  + V)  t 


(35) 


The  quantity  A0  is  the  concentration  difference  at  the  onset  of  a linear  concentration 
profile.  Based  on  Eqn.  35,  a semi-log  plot  of  the  concentration  difference  between  the 
two  vessels  should  appear  as  a straight  line  under  ideal  conditions;  the  magnitude  of 
the  slope  being  proportional  to  the  apparent  diffusion  coefficient  Da.  Deviations  from  a 
straight  line  will  indicate  behavior  that  cannot  be  modeled  by  Fick’s  law  with  a constant 
apparent  diffusion  coefficient. 


4.2.4  Formation  Factor  Measurement 

The  formation  factor  measurements  were  performed  using  KC1  as  the  pore  solution  be- 
cause the  conductivity  of  a few  standard  solutions  are  known  to  a high  precision  [37].  A 
range  of  concentrations  was  used  to  ensure  that  the  surface  conduction  component  was 
properly  accounted  for.  For  the  conductivity  measurement,  the  specimen  was  vacuum 
saturated  with  potassium  chloride  solution  and  then  mounted,  using  rubber  o-rings  and 
clamps,  between  two  glass  vessels,  each  containing  a platinum  electrode.  The  setup  is 
shown  schematically  in  Fig.  10(a).  The  apparatus,  in  the  absence  of  a specimen  and 
holder,  had  a conductivity  cell  constant  of  0.3567  cm-1. 

To  determine  the  formation  factor  of  the  saturated  specimen,  the  cell  was  filled  with 
the  same  solution  as  the  saturation  solution  and  then  allowed  to  thermally  equilibrate  in 
the  environmental  chamber.  The  direct  current  (dc)  resistance  of  the  specimen  and  cell 
was  determined  using  a commercial  impedance  spectrometer  that  sampled  frequencies 
between  10  Hz  and  1 MHz  [22].  The  bulk  conductivity  of  the  specimen  was  calculated 
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(a)  (b) 

Figure  10:  The  conductivity  cell  apparatus(a)  and  the  diffusion  divided  cell  apparatus 
(b).  The  schematics  depict  the  configuration  of  the  two  cylindrical  glass  vessels  on  either 
side  of  a mounted  specimen.  The  system  is  sealed  using  rubber  o- rings;  the  clamps  are  not 
shown.  The  apparatus  differ  only  in  the  vertical  platinum  electrodes  in  the  conductivity 
cell.  The  diameter  of  the  specimen,  the  glass  vessels,  and  the  platinum  electrodes  are 
approximately  50  mm. 

from  the  cell  constant  and  the  specimen  geometry.  The  formation  factor  was  calculated 
from  the  ratio  of  the  conductivity  of  the  KC1  solution  to  the  calculated  bulk  conductivity. 

Pore  solution  concentrations  of  0.01  mol-kg-1,  0.10  mol-kg-1,  and  1.00  mol-kg-1  were 
used  to  assess  the  contribution  from  surface  conduction.  Due  to  a small  surface  conduc- 
tion contribution,  the  formation  factor  increased  with  increasing  pore  solution  conductiv- 
ity, converging  to  a fixed  value  with  increasing  concentration.  The  specimen  conductivity 
measured  using  the  0.01  mol-kg-1  solution  was  approximately  85  % of  the  value  using  the 
1.0  mol-kg-1  solution,  and  the  specimen  conductivity  using  the  0.10  mol-kg-1  solution 
was  approximately  98  % of  the  value  using  the  1.0  mol-kg-1  solution.  Therefore,  the 
formation  factor  calculated  at  1.0  mol-kg-1  was  used  as  the  best  estimate.  Since  the 
change  in  formation  factor  was  only  2 % for  a ten  fold  increase  in  concentration,  the  true 
value  is  probably  not  more  than  a fraction  of  percent  larger  than  the  values  reported. 

Since  the  ceramic  specimens  are  the  result  of  a controlled  commercial  process,  they 
have  relatively  little  specimen  to  specimen  variation.  The  four  specimens  used  in  this  ex- 
periment had  formation  factors  ranging  from  10.6  to  10.9.  These  values  were  determined 
after  each  specimen  was  allowed  sufficient  time  to  reach  thermal  equilibrium,  and  then 
the  values  of  the  dc  resistance  varied  by  less  than  1 %. 

4.3  Formation  Factor  vs.  Diffusion  Coefficient 

Upon  completion  of  the  conductivity  measurement  to  determine  the  formation  factor, 
the  specimen  was  then  saturated  with  KI  and  placed  into  a divided  cell  apparatus  shown 
schematically  in  Fig.  10(b).  The  objective  was  to  determine  the  value  for  the  formation 
factor  F required  in  Eqn.  32  in  order  for  the  calculated  iodide  concentrations  to  agree 
with  the  experimental  values. 

The  measured  values  of  A for  the  iodide  concentration  are  shown  as  symbols  in  Fig.  11. 
For  each  of  the  four  systems,  the  concentration  of  the  KI  is  equal  to  the  concentrations 
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Figure  11:  Concentration  difference  across  each  specimen  as  a function  of  time  normalized 
by  the  vessel  volumes.  The  experimental  values  are  shown  as  Filed  symbols,  the  estimates 
by  4SIGHT  are  shown  as  solid  curves.  The  measurement  uncertainties  would  appear  as 
the  same  size  as  the  symbols,  so  are  omitted  for  visual  clarity.  The  value  of  A for  the  1.0 
mol/L  KC1  system  is  divided  by  ten  in  order  to  appear  on  the  same  scale  as  the  other 
data. 

shown  in  the  figure  key.  The  estimated  measurement  uncertainties  are  approximately 
the  size  of  the  symbols,  and  are  not  shown  as  error  bars  for  reasons  of  visual  clarity. 
The  KC1  and  the  NaCl  systems  behaved  similarly;  the  values  of  A for  the  1.0  mol/L 
KC1  system  were  divided  by  ten  so  that  they  could  be  included  on  the  same  plot.  The 
1.0  mol/L  KC1  solution  has  a similar  behavior  to  the  0.1  mol/L  KC1  system  because 
the  self  diffusion  coefficients  Ds  for  K+,  Cl~,  and  I-  are  nearly  equal  to  one  another. 
The  estimated  diffusion  potential  t/’d  (Eqn.  32  across  the  entire  specimen  for  both  of  the 
KC1  systems  was  less  than  1 mV,  and  was  less  than  5 mV  for  the  NaCl  system.  The 
KOH  system  showed  a marked  difference  in  behavior.  One  reason  for  this  is  that  the 
self  diffusion  coefficient  Ds  for  OH~  is  significantly  greater  than  that  of  the  other  ions 
present,  resulting  in  a calculated  diffusion  potential  of  approximately  16  mV.  A diffusion 
potential  of  this  magnitude  could  have  a noticeable  effect  on  the  transport  of  the  iodide 
ions. 

The  calculated  slopes  of  the  experimental  data  from  semi-logarithmic  plots  are  shown 
in  Table  4.  The  estimated  standard  deviations  shown  are  typically  less  than  3 % of  the 
slope  value,  suggesting  that  the  constant  gradient  assumption  is  valid  for  these  systems. 
In  fact,  the  systems  reached  a constant  gradient  state  at  a relatively  early  age.  Analysis 
using  the  4SIGHT  computer  program  suggests  that  the  constant  gradient  condition  is 
achieved  in  less  than  12  h. 

Also  shown  in  Table  4 are  the  values  for  the  ratio  Ds /Dt,,  using  the  iodide  value  for 
Ds  (2.045  x 10~5  cm2  s'1  [23]).  This  ratio  represents  an  incorrect  use  of  the  formation 
factor  to  determine  the  apparent  bulk  diffusion  coefficient.  Since  the  quantity  Ds  in  this 
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Table  4:  The  values  for  the  slope  ( ADb/L ) of  the  experimental  data  shown  in  Fig.  11. 
Also  shown  is  the  ratio  Ds/Db , using  Ds  for  iodide.  The  uncertainties  shown  for  the 
slopes  are  the  estimated  standard  deviation  reported  by  the  statistical  software,  and  also 
characterize  the  uncertainty  in  the  ratio  Ds / Db  reported. 


System 

ADb/L 

(cmOr1) 

Ds/Db 

KC1 

- 0.1  mol/L 

0.2004  ±0.0032 

11.1 

NaCl 

- 0.1  mol/L 

0.21 18±0. 0032 

10.3 

KOH 

- 0.1  mol/L 

0.2381  ±0.0047 

9.3 

KC1 

1 .0  mol/L 

0.2079±0.0050 

10.7 

ratio  is  a constant,  the  ratio  is  inversely  proportional  to  the  apparent  bulk  diffusivity  Dh . 
This  ratio  does  not  reflect  the  actual  formation  factor,  however,  because  the  iodide  self 
diffusion  coefficient  within  the  pore  solutions  is  not  equal  to  Ds . Also,  it  is  clear  that 
future  changes  in  the  pore  solution  will  lead  to  changes  in  the  apparent  bulk  diffusion 
coefficient  of  iodide  in  these  systems,  while  the  formation  factor  is  nearly  equal  for  all 
systems. 

The  values  of  the  formation  factor  F were  determined  from  the  experimental  data  us- 
ing the  4SIGHT  implementation  of  Eqn.  32,  and  are  shown  in  Table  5.  labelled  Fstm.  Also 
shown  in  the  table  are  the  values  of  the  formation  factor  calculated  from  the  impedance 
spectroscopy  measurements,  labelled  FIS.  The  uncertainty  in  Fts  reflects  the  variation 
in  the  dc  resistance  measurement  as  mentioned  previously.  To  estimate  F in  Eqn.  32, 
its  value  was  adjusted  until  the  calculated  values  of  A agreed  with  the  measured  values; 
this  value  is  referred  to  here  as  Fsim.  The  values  of  Fsim  shown  in  Table  5 were  used  to 
calculate  values  for  A.  and  these  values  of  A are  plotted  in  Fig.  11,  denoted  by  the  lines. 
The  experimental  data  and  the  calculated  values  for  the  KC1  and  the  NaCl  systems  were 
all  nearly  linear,  suggesting  that  the  linear  Fick’s  law  for  diffusive  transport  is  valid  for 
these  systems.  The  values  of  A for  the  KOH  system  are  easily  distinguished  from  the 
other  systems. 

The  calculated  values  of  Fsim  shown  in  Table  5 were  consistant  with  the  measured 
values  FIS ■ The  values  of  Fsim  varied  by  approximately  7 %,  compared  with  the  18  % 
variation  in  the  values  of  Ds / Db.  The  differences  between  the  values  of  Fsim  and  Fjs 
were  less  than  3 % for  the  KC1  and  the  NaCl  system,  and  less  than  8 % for  the  KOH 
system.  The  consistant  bias  of  calculated  values  of  Fsim  being  larger  than  Fls  may  be 
attributed  to  the  method  of  using  a single  measurement  for  the  specimen  conductivity, 
rather  than  using  a more  rigorous  approach. 

Sighting  along  the  KOH  data  in  Fig.  11  reveals  noticeable  curvature.  It  is  interesting 
to  note  that  the  computed  output  (solid  curve)  also  exhibits  this  nonlinear  behavior.  This 
suggests  that  this  nonlinear  behavior  is  due  to  effects  of  the  pore  solution  chemistry  since 
output  from  4SIGHT  indicates  that  the  iodide  concentration  profile  across  the  sample  is 
stable  within  12  h.  The  calculated  concentration  profile  of  iodide,  however,  is  not  linear 
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due  to  the  diffusion  potential 


Table  5:  Measured  and  calculated  formation  factors  from  impedance  spectroscopy  (T75), 
computer  simulation  ( Fsirn ) , and  apparent  diffusivity  ( Ds / Db ) using  Ds  for  iodide. 


System 

Fis 

F 

1 sun 

D‘/Db 

KC1  - 0.1  mol/L 

10.7T0.2 

10.9 

11.1 

Nad  - 0.1  mol/L 

10.9T0.2 

11.2 

10.5 

KOI I - 0.1  mol/L 

10.6T0.2 

11.4 

9.3 

KC1  - 1.0  mol/L 

10.7T0.2 

10.6 

10.7 

This  electro-chemical  effect  of  the  KOH  system  is  revealed  in  Table  5.  As  shown  in 
Table  4,  the  apparent  diffusion  coefficient  of  iodide  in  this  system  is  considerably  greater 
than  that  for  the  other  systems.  The  computer  calculation  reveals  that  the  required 
formation  factor  Fsim  is  comparable  to  the  values  for  the  other  test  solutions.  This 
demonstrates  the  effect  of  using  the  apparent  diffusion  coefficient  Db  to  characterize  a 
microstructure.  For  KOH  solution,  the  apparent  diffusion  coefficient  describes  how  the 
iodide  ion  behaves  in  the  presence  of  KOH,  but  does  not  characterize  its  behavior  in  the 
presence  of  other  test  solutions.  Similarly,  it  does  not  necessarily  characterize  how  other 
ions  behave  in  the  same,  or  similar,  microstructure. 

Since  the  pore  solution  of  cementitious  systems  is  typically  alkaline,  the  results  for  the 
KOH  system  have  direct  relevance  to  the  prediction  of  ion  transport  in  portland  cement 
systems.  The  pore  solution  ionic  strength  in  cementitious  systems  can  be  nearly  ten  times 
greater  than  the  0.1  mol/L  KOH  system  studied  here.  Further,  there  will  be  number  of 
additional  ions  present,  with  a corresponding  number  of  self  diffusion  coefficients.  This 
raises  the  question  of  the  correct  method  for  characterizing  the  microstructure  of  these 
systems.  Either  estimating  the  formation  factor  from  diffusion  data  or  estimating  the 
diffusion  coefficient  from  a formation  factor  measurement  (conduction  test)  will  require 
knowledge  of  the  pore  solution  chemistry.  For  a formation  factor  measurement  that  uses 
pore  solution  extraction,  a chemical  analysis  of  the  extracted  pore  fluid  would  be  a logical 
extension  of  the  measurement  procedure.  However,  diffusion  measurements  typically  do 
not  reveal  the  chemical  makeup  of  the  pore  solution,  and  so  further  experimentation  may 
be  required  for  experimental  programs  based  upon  diffusion  measurements  to  characterize 
the  microstructure  of  cementitious  systems. 

4,3.1  Multi-Species 

The  results  for  the  multi-species  experiments  are  divided  between  two  series.  I11  the 
first  series,  the  kinds  of  ionic  species  present  were  kept  constant  but  the  initial  iodide 
concentration  differences  were  varied.  In  the  second,  the  initial  iodide  concentration 
were  held  constant,  but  the  kinds  of  other  species  present  were  varied.  The  first  series 
attempts  to  demonstrate  the  effect  of  concentration.  The  second  attempts  to  demonstrate 
the  effect  of  speciation. 
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Figure  12:  Concentration  difference  A of  iodide  (filled  symbols)  for  the  KC1/KI  and  the 
NaCl/KI  systems;  the  time  t is  scaled  by  the  vessel  volumes  v\  and  u2.  The  quantity  A0 
is  the  initial  KI  concentration.  The  lines  shown  are  the  result  of  the  linear  regression 
analyses  reported  in  Table  6.  The  measurement  uncertainties  are  approximately  the  size 
of  the  symbols  used. 


The  experiments  with  varying  concentration  were  performed  with  either  potassium 
chloride  or  sodium  chloride  as  the  counter-diffusing  electrolyte,  with  respect  to  potasium 
iodide.  Plots  of  the  experimental  results  are  shown  in  Figs.  12.  For  both  the  KC1/KI  and 
the  NaCl/KI  systems,  the  data  exhibit  virtually  linear  behavior  on  the  semi-logarithmic 
plots.  The  apparent  linear  behavior  suggests  that  these  systems  behave  in  a nearly 
“ideal”  manner,  and  can  be  characterized  accurately  by  Fick’s  law  in  Eqn.  10.  Due  to 
the  nearly  equal  self-diffusion  coefficients  for  potassium,  chloride,  and  iodide,  there  is 
only  a relatively  small  electrical  diffusion  potential  across  the  sample. 

The  magnitude  of  the  slope  calculated  from  the  data  at  each  concentration  is  reported 
in  Table  6.  The  uncertainties  reported  are  the  estimated  standard  deviations  reported 
by  the  regression  software,  and  have  the  corresponding  coverage  factor  [40]  k = 1.  For 
the  three  concentrations,  the  coefficient  of  variation  of  the  estimated  slopes  was  5 % for 
the  KC1/KI  systems,  and  10  % for  the  NaCl/KI  systems.  While  the  differences  among 
slopes,  being  greater  than  the  estimated  uncertainties,  suggest  that  these  variations  are 
statistically  significant,  for  practical  use,  these  variations  are  minor.  Any  chemical  reac- 
tions that  might  occur  within  the  cementitious  matrix  can  alter  the  diffusion  coefficient, 
easily  overwhelming  a 20  % uncertainty  in  the  apparent  diffusion  coefficient. 

The  experiments  with  constant  KI  concentration  and  varying  speciation  were  per- 
formed with  KC1,  NaCl,  KOH,  K2C03,  CaCl2,  H2S04,  and  MgSCK  as  the  counter- 
diffusing  electrolyte.  In  all  cases,  the  initial  concentration  of  KI  was  0.10  mol/L.  The 
concentration  of  the  counter-diffusion  electrolyte  was  chosen  to  minimize  an  osmotic  pres- 
sure gradient  across  the  sample;  it  is  not  asserted  here  that  the  values  used  eliminated 
an  osmotic  pressure  gradient. 
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Table  6:  Absolute  value  of  the  slopes  calculated  from  the  data  for  the  KC1/KI  and 
the  NaCl/KI  systems  in  Figs.  12,  regressed  to  a linear  model.  The  uncertainties  shown 
are  the  estimated  from  standard  deviations  reported  by  the  regression  software,  with  a 
corresponding  coverage  factor  of  k = 1. 


Species 

concentration 

(mol/L) 

| Slope  | 
(cm3/h) 

KC1/KI 

0.010/0.010 

0.100/0.100 

1.000/1.000 

0.1923T0.0039 

0.1978T0.0029 

0.2101T0.0024 

NaCl/KI 

0.010/0.010 

0.100/0.100 

1.024/1.000 

0.2365T0.0051 

0.2120T0.0024 

0.1931T0.0059 

The  experimental  data  are  shown  in  Figs.  13,  and  the  analyses  of  the  data  are  divided 
between  the  “short”  time  and  the  “long”  time  results.  In  Fig.  13(a),  a linear  model  was 
regressed  to  the  data  for  which  t(v f1  + v2  l)  < 8 h/cm3.  The  results  from  the  linear 
regression  analyses  are  shown  in  Table  7.  The  uncertainties  reported  are  the  estimated 
standard  deviations  reported  by  the  regression  software,  and  have  the  corresponding  cov- 
erage factor  [40]  k = 1.  Based  on  these  results,  the  CaCC/KI  system  and  the  K0CO3/KI 
system  represent  the  extrema,  and  are  shown  in  Fig.  13(b)  for  all  time,  along  with  the 
KC1/KI  system,  which  is  an  example  of  nearly  ideal  behavior,  as  a means  to  constrast 
the  results. 

Table  7:  Magnitude  of  the  slope  for  data  shown  in  Fig.  13(a),  regressed  to  Eqn.  35 
for  t(i T1  + n^1)  < 8 h/cm3.  For  each  system,  the  I\I  concentration  was  0.10  rnol/L. 
Also  shown  are  the  values  for  the  formation  factor  F used  to  produce  the  solid  curves 
Fig.  13(b).  The  uncertainties  shown  are  the  estimated  from  standard  deviations  reported 
bv  the  regression  software,  with  a corresponding  coverage  factor  of  h — 1. 


Counter  Solution 

concentration 

(mol/L) 

|Slope| 

(cm3/h) 

F 

k2co3 

0.0769 

0.1581T0.0014 

11.7 

H2SO4 

0.0859 

0. 1688T0.0075 

9.6 

MgSCC 

0.1660 

0.1867T0.0022 

13.5 

KC1 

0.1000 

0.1978T0.0029 

11.4 

Nad 

0.1000 

0.2120T0.0024 

11.8 

KOH 

0.0879 

0.2333T0.0034 

11.2 

CaCl2 

0.0731 

0.2801  ±0.0089 

12.8 
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Figure  13:  Measured  concentration  difference  A of  iodide  (filled  symbols)  for  constant 
0.10  mol/L  KI  concentration  and  varying  counter-diffusing  species;  the  time  f is  scaled 
by  the  vessel  volumes  V\  and  V2 . In  (a)  results  from  linear  regression  (solid  lines)  for 
t,(v f1  + U71)  < 8 h/cm3;  the  calculated  slopes  are  reported  in  Table  7.  In  (b)  the  electro- 
diffusion equation  is  used  to  estimate  the  formation  factor  F,  the  values  of  which  are 
shown  in  Table  7.  The  error  bars  represent  an  estimated  standard  deviation. 


The  computed  slopes  given  in  Table  7 indicate  the  effect  of  speciation  on  the  observed 
diffusion  coefficient.  The  magnitude  of  the  slopes  shown  in  the  table  varies  by  nearly  a 
factor  two.  This  contrasts  with  the  previous  20  % variation  in  slopes  due  to  concentration. 

The  CaClo  and  the  K2CO3  data  are  shown  in  Fig.  13(b)  to  delineate  the  envelop  over 
which  the  concentrations  could  vary  over  long  times.  Note  that  for  long  times  the  slopes 
of  the  curves  are  not  constant.  While  the  change  in  slope  for  the  K2CO3/KI  system  is 
gradual,  the  change  in  slope  for  the  CaCF/KI  system  is  dramatic  by  comparison.  There- 
fore, it  is  inappropriate  to  characterizing  iodide  transport  in  the  CaCl2/KI  system  using 
Fick’s  law  and  a constant  apparent  diffusion  coefficient.  Further,  the  apparent  iodide  dif- 
fusion coefficient  decreases  slightly  in  the  K2CO3/KI  system  and  increases  greatly  in  the 
CaCF/KI  system.  Any  “adjustments”  to  an  apparent  diffusion  coefficient  to  account  for 
one  system  would  certainly  not  apply  to  another.  This  is  analogous  to  determining  the 
apparent  diffusion  coefficient  under  one  set  of  conditions  and  using  that  to  extrapolate 
to  behavior  in  the  held  with  varying  species  exposure. 

By  contrast,  the  data  for  each  system  were  analyzed  using  the  4SIGHT  implementation 
of  Eqn.  32.  The  porosity  remained  fixed  at  0.26  and  the  formation  factor  remained  the 
only  free  parameter.  The  value  of  the  formation  factor  F required  to  best  approximate 
the  measured  data,  for  all  times,  is  shown  in  Table  7 for  each  system.  The  coefficient  of 
variation  for  these  values  was  approximately  10  %.  The  corresponding  results  are  shown 
as  solid  curves  111  Fig.  13(b).  The  curves  for  the  other  systems  are  not  shown  for  clarity. 
Further,  remaining  systems  had  only  minor  curvature  by  comparison. 

The  important  feature  of  the  electro-diffusion  equation  is  that  it  can  capture  the 
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Figure  14:  The  ratio  of  the  apparent  diffusion  coefficient  Da  to  the  microstructural 
diffusion  coefficient  Du  (as  calculated  by  4SIGffT  using  Eqn.  14)  of  iodide  for  the  various 
counter  diffusing  electrolytes,  as  a function  of  time  t. 


nonideal  behavior  of  these  systems.  The  apparent  diffusion  coefficient  can  be  calculated 
at  any  time  from  the  slope  of  the  curve  based  on  the  electro-diffusion  equation.  This 
calculated  apparent  diffusion  coefficient  can  then  be  compared  with  the  microstructural 
diffusion  coefficient.  The  results  from  this  comparison  are  shown  in  Fig.  14  for  all  the 
systems  considered  in  Fig.  14. 

The  curves  shown  in  Fig.  14  express  quantitatively  the  deficiency  in  characterizing 
multi-species  diffusion  with  a single  apparent  diffusion  coefficient.  Accurately  charac- 
terizing the  diffusive  transport  of  iodide  in  these  nonreactive  systems  using  Fick’s  law 
requires  that  the  apparent  diffusion  coefficient  vary  in  time.  For  the  CaCl2/KI  system, 
the  apparent  diffusion  coefficient  has  increased  by  a factor  of  5 at  60  d. 

The  surprising  result  was  the  behavior  of  the  CaC^/Kl  system  at  very  long  time. 
The  curve  for  the  CaCl2/Kl  system  in  Fig.  13(b)  appears  to  be  diverging  at  long  times. 
Analysis  of  the  experimental  setup  using  the  electro-diffusion  equation  predicts  that  the 
concentration  difference  between  the  vessels  will  become  negative ; this  is  equivalent  to 
a negative  apparent  diffusion  coefficient  as  the  flux  of  iodide  ions  is  in  the  direction  of 
the  low  concentration  to  high  concentration.  The  concentration  differences  A for  the 
CaC^/KI  system  are  shown  in  Fig.  15(a)  on  a linear  scale.  After  approximately  70  d, 
the  experimentally  measured  concentration  difference  was,  in  fact,  negative.  As  a means 
of  comparison,  the  corresponding  data  for  the  K2CO3/KI  system  are  shown  in  Fig.  15(b). 

Also  shown  in  Figs.  15(a)  and  15(b)  are  the  relative  concentration  differences  for 
all  the  species  present,  as  calculated  by  Eqn.  35.  For  the  CaCl2/KI  system,  the  “slow” 
species  is  the  Ca2+  ion.  For  the  K2CO3/KI  system,  it  is  the  CO2-  ion.  The  effect  of  these 
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Figure  15:  Experimental  concentration  difference  A of  each  diffusing  species  in  the  (a) 
CaCl2/KI,  and  the  (b)  K2C03/KI  systems,  each  as  a function  of  time  t.  The  experimental 
values  are  shown  as  filled  symbols,  and  the  curves  are  the  values  calculated  by  4SIGHT. 
The  error  bars  represent  an  estimated  standard  deviation. 


two  ions  is  to  greatly  control  the  macroscopic  diffusion  potential  across  the  specimen.  The 
presence  of  the  Ca2+  generates  a diffusion  potential  that  is  positive  with  respect  to  the 
KI  side.  Correspondingly,  the  C03~  ion  generates  a negative  potential  with  respect  to 
the  KI  side.  The  result  is  a macroscopic  diffusion  potential  that  accelerates  the  iodide 
ion  in  the  CaCl2  system,  and  retards  it  in  the  K2C03  system.  In  the  CaCl2  system, 
the  macroscopic  diffusion  potential  is  sufficient  to  drive  iodide  from  the  KI  vessel  to  the 
CaCl2  vessel,  even  though  the  iodide  concentration  in  the  KI  vessel  is  less  than  that  of  the 
CaCl2  vessel.  In  time,  the  diffusion  potential  decreases  and  the  concentration  difference 
A approaches  zero. 

4.4  Summary 

Since  each  of  the  degradation  mechanisms  considered  by  4SIGHT  are  controlled  by  the 
transport  of  ionic  species  throughout  the  concrete,  it  is  imparitive  that  the  transport 
portion  of  the  4SIGHT  computer  program  be  as  accurate  as  possible.  The  first  step  was  to 
validate  the  transport  equation  Eqn.  32  used  by  4SIGHT.  It  was  shown  that  the  formation 
factor  used  in  Eqn.  32  is  equivalent  to  the  formation  factor  that  can  be  calculated  from 
conductivity  measurements.  The  second  step  was  to  demonstrate  that  the  transport 
equation  characterizes  the  transport  of  all  the  ions  present  in  the  electrolyte,  and  that 
the  only  two  parameters  required  are  the  porosity  and  formation  factor.  By  characterizing 
the  transport  of  all  the  species  present,  and  the  resulting  macroscopic  diffusion  potential, 
the  4SIGHT  computer  program  successfully  predicted  the  occurance  of  an  experimental 
setup  that  could  exhibit  a negative  apparent  diffusion  coefficient  (based  on  Fields  law). 
Although  only  iodide  concentration  measurements  were  performed,  the  wide  variation 
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of  behavior  of  the  different  experimental  setups  suggests  that  the  complete  system  was 
modeled  accurately.  Future  experiments  where  all  diffusion  species  concentrations  are 
measured  are  needed  to  confirm  this. 


5 ' Parameter  Uncertainty 

As  the  complexity  of  service  life  estimation  increases,  it  is  clear  that  deterministic  cal- 
culations have  diminishing  utility.  One  cannot  say  with  absolute  certainty,  what  the 
performance  of  a facility  will  be  over  a period  of  decades  or  centuries.  Rather,  one  can 
only  offer  probabilistic  statements  of  possible  outcomes. 

To  facilitate  a probabilistic  risk  assessment,  4SIGHT  allows  the  user  to  specify  input 
parameters  with  both  a mean  and  a level  of  uncertainty.  These  uncertainty  parameters 
are  then  used  to  generate  random  deviates  based  on  a random  number  generator.  Multi- 
ple service  life  calculations  are  performed,  each  using  a new  random  deviate.  The  failure 
times  for  each  calculation  is  recorded.  The  results  are  then  reported  as  a cumulative 
distribution  of  failure  times. 

5.1  Random  Deviates 

The  basis  for  the  random  deviate  calculation  lies  in  the  ability  to  generate  a random 
deviate  from  the  uniform  distribution  over  the  interval  (0,1).  This  is  accomplished  using 
a random  number  generator  (see  Ch.  7 of  Ref.  [41]).  The  resulting  uniform  random 
deviate  i]}  is  then  used  to  calculate  a random  deviate  Xt  from  an  arbitrary  probability 
density  function  f(x): 

f f(x)dx  = i],  (36) 

J — OG 

The  random  number  generator  used  by  4SIGHT  is  from  Knuth  [42],  and  referred  to  as 
ran3  in  Numerical  Recipes  [41].  It  has  been  shown  that  this  random  number  generator 
is  sufficiently  robust  for  calculations  involving  1 0 ' deviates. 

5.2  Distributions 

Four  probability  density  functions  f(x)  are  provided  to  the  user  by  4SIGHT.  These  four 
distributions  can  be  used  to  characterize  random  deviates  over  a broad  range  of  uncer- 
tainty. 

5.2.1  Delta  Function 

The  delta  function  (see  Refs.  [43,44])  is  equivalent  to  having  a deterministic  value.  The 
delta  function  is,  in  a practical  sense,  a Gaussian  distribution  with  zero  variance.  The 
symbol  typically  used  for  the  delta  function  is  5{x  — x0),  and  a plot  of  the  function  would 
consist  of  a vertical  line  at  i = xQ.  A delta  function  has  the  following  mathematical 
definitions  [44]: 


/ S(x)  dx  = 1 

J — OO 

g{x 0)  a < x0  < b 


g(x)S(x  — x0)  dx  = 


0 


otherwise 


(37) 


The  probability  density  function  f(x)  is  written  mathematically  as  follows: 

f(x)  = 5(x  — .t0)  — oo  < x < Too  (38) 

Every  random  deviate  xl  will  have  the  same  value  t0: 

x i = t0 


5.2.2  Uniform  Distribution 

The  uniform  distribution  allows  one  to  specify  a range  over  which  values  can  vary;  values 
cannot  be  used  that  lie  outside  the  specified  width  of  the  distribution.  As  its  name 
implies,  random  deviates  are  chosen  uniformly  over  the  specified  range,  characterized 
by  a center  at  t0  and  a width  w.  Values  are  chosen  uniformly  between  xQ  — w/2  and 
x0  + w/2 : 

„ / , 1 w w , 

/( x)  = — To-  — <T<T0  + — (39) 

w 2 2 

Using  the  uniform  random  deviate  r\t  over  the  interval  (0,1),  the  random  deviate  xt 
is  calculated  from  the  center  x0  and  width  w of  the  distribution: 

Xi  = x0  + w (ru  - - \ 


5.2.3  Gaussian  Distribution 

The  Gaussian  distribution  is  familiar.  For  a mean  /i  = xQ  and  variance  ex2,  the  probability 
density  function  can  be  stated  mathematically  as  follows: 


f(x) 


(40) 


Using  standard  software  routines  [41],  a random  deviate  A is  calculated  from  the  Gaussian 
distribution  with  zero  mean  and  unit  standard  deviation.  The  normal  random  deviate  G 
can  then  be  used  to  calculate  a random  deviate  from  a Gaussian  distribution  with  mean 
xQ  and  variance  a2: 

Xi  = T0  + cr<G 


5.2.4  Lognormal  Distribution 

An  advantage  of  the  lognormal  distribution  is  that  it  can  characterize  uncertainty  over 
orders  of  magnitude.  Although  an  entire  family  of  n-th  order  logarithmic  distributions 
exist  [45],  the  basic  feature  of  a lognormal  distribution  is  that  the  logarithms  of  the  devi- 
ates are  normally  distributed.  This  fact  can  be  exploited  to  simplify  deviate  generation. 

Consider  experimentally  determined  quantities  yx  that  are  the  result  of  a measurement 
process,  and  can  be  characterized  by  a lognormal  distribution  (e.g.,  hydraulic  conductiv- 
ity). Let  the  variable  z%  represent  the  logarithms  (base  10)  of  these  values:  zt  — log10(y;). 
Let  x0  denote  the  average  of  the  values  of  zt.  Let  a represent  the  standard  deviation  of 
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the  values  zt.  Given  an  normal  random  deviate  G with  zero  mean  and  unit  variance,  a 
logarithmically  distributed  random  deviate  xt  can  be  generated  in  the  following  way: 


(41) 


5.3  Output:  Cumulative  Distribution 

For  the  Monte  Carlo  calculations,  the  user  will  specify  a number  N of  iterations.  For  each 
iteration,  new  random  deviates  will  be  calculated  and  the  time  at  which  failure  occurs 
will  be  calculated;  in  the  absence  of  failure,  the  program  returns  the  time  limit  set  by  the 
user.  Let  each  failure  time  be  denoted  by  6,  where  l<  i < N . The  probability  density 
function  g(t. ),  giving  the  probability  that  the  service  life  of  a computational  iteration 
terminated  between  t and  t + dt,  can  be  expressed  as  a sum  of  delta  functions: 


(42) 


The  cumulative  probability  function  G(t)  is  the  integral  of  the  probability  density 
function  g(t): 


(43) 


Therefore,  under  the  parameter  uncertainty  conditions  specified  by  the  user,  the  proba- 
bility that  failure  will  occur  before  time  r is  simply  G(r). 

This  is  a mathematically  rigorous  means  of  expressing  a straightforward  concept.  An 
example  will  clarify  the  idea  further:  Consider  a service  life  calculation  performed  with 
N = 10  iterations,  with  each  iteration  using  different  input  parameters,  resulting  in  the 
following  service  life  termination  times: 


1  120 

2 324 

3 80 

4 260 

5 250 

6 210 

7 500 

8 170 

9 500 

10  280 


The  corresponding  cumulative  distribution  function  G(t)  is  given  in  Fig.  16.  These  data 
can  then  be  used  within  a larger  service  life  calculation  composed  of  a number  of  elements, 
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Figure  16:  The  cumulative  distribution  function  G(t). 


of  which  this  is  only  one.  A random  failure  deviate  can  be  generated  from  these  data 
by  picking  a uniform  random  deviate  r/2  over  the  interval  (0,1)  and  solving  for  t,r  from 
G{U)  = T)i . 
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6 Significance  of  Cracks 

Most  degradation  processes  are  due  to  the  presence  of  some  specific  ion  within  the  cement 
paste  microstructure.  These  deleterious  ions  typically  penetrate  the  microstructure  by 
diffusion  and/or  advection.  To  obtain  an  understanding  of  the  relative  magnitude  of 
these  mechanisms,  consider  their  contribution  to  the  ionic  flux.  For  diffusivity  D , ion 
concentration  c within  the  capillary  pores,  and  pore  solution  velocity  up,  the  flux  j is  a 
sum  of  contributions  from  diffusion  and  advection: 

j = — D'Ve  + cup 

V n 

= -DVc  + c— 

<t> 

= — DVc+c — VP  (44) 

4>k 

Here,  the  permeability  k is  the  bulk  permeability  that  is  often  reported.  The  second  line 
in  Eqn.  44  represents  the  pore  velocity  as  a function  of  the  Darcy  (bulk)  velocity  vD 
and  capillary  porosity  0.  The  third  line  contains  the  substitution  for  the  Darcy  velocity 
as  a function  of  the  permeability  k , pore  fluid  viscosity  //,  and  hydraulic  pressure  P. 
For  a 0.45  waterxement  ratio  (w/c)  concrete,  the  relevant  material  parameters  can  be 
estimated  from  the  previous  4SIGHT  report  [1]  and  are  shown  in  Table  8. 

Table  8:  Material  parameters  a concrete  with  w/c  = 0.45,  estimated  from  equations  in 
Ref.  [1].  (Default  degree  of  hydration  a — 0.85) 


Parameter 

Symbol 

Value 

diffusivity 

D 

7x10“ 12  m2s_1 

permeability 

k 

lx  10“ 19  m2 

viscosity 

F 

10~3  Pa  s 

porosity 

4> 

0.11 

The  gradients  in  Eqn.  44  can  be  estimated  from  “order  of  magnitude”  values.  The 
length  scale  is  approximately  1 m,  the  concentrations  on  the  order  of  1 mol,  and  the 
hydraulic  pressure  a tenth  of  an  atmosphere  (104  Pa).  The  contributions  to  the  flux  are 
7xl0-12  for  the  diffusive  transport  and  9xl0-12  for  the  advective  transport.  Therefore, 
this  scenario  represents  an  approximate  equality  between  the  relative  contributions  to 
the  overall  flux.  Changes  in  either  the  concentration  or  the  hydraulic  gradients  would 
change  the  contributions,  respectively.  Since  cracks  in  concrete  have  permeabilities  that 
are  orders  of  magnitude  greater  than  the  permeability  of  uncracked  concrete,  their  impact 
on  the  ionic  flux  can  be  dramatic. 

6.1  Internal  Body  Force 

The  previous  report  [1]  contained  a model  for  advective  transport  through  uncracked 
concrete  that  included  a term  for  the  body  force  of  pg  ( p is  the  fluid  density  and  g is 
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Figure  17:  Schematic  of  water  sorption  in  concrete.  The  hashed  region  represents  a cross 
section,  with  pore  radius  R and  liquid-air  surface  tension  7,  with  its  bottom  surface 
under  water.  To  the  right  is  a plot  of  the  hydraulic  pressure  P as  one  traverses  vertically 
through  the  pore.  Negative  pressure  indicates  water  in  tension. 


the  gravitational  acceleration),  adding  a constant  to  the  pressure  gradient.  This  term 
contributed  to  the  overall  advect.ive  transport  through  the  concrete.  While  this  addition 
is  accurate  for  a large  saturated  specimen  with  water  at  both  ends,  it  is  not  accurate 
when  either  end  is  exposed  to  an  unsaturated  environment.  Exposure  to  an  unsaturated 
environment  will  lead  to  the  formation  of  menisci  in  the  cement  paste  pore  system. 
Under  many  possible  conditions,  the  pressures  developed  by  these  menisci  are  sufficient 
to  prevent  the  water  from  “draining”  out  of  the  concrete  due  to  gravity. 

6.2  Sorptivity 

Typical  capillary  sorptivity  tests  demonstrate  that  concrete  is  not  a self-draining  material. 
A dry  concrete  specimen  placed  onto  a free  surface  of  water  will  draw  the  water  into  the 
concrete  due  to  capillary  sorption.  If  allowed  to  contact  the  water  for  a sufficiently 
long  period  of  time,  and  assuming  the  height  of  the  concrete  specimen  is  less  than  10 
m (hydraulic  head  of  one  atmosphere),  the  concrete  will,  given  sufficient  time,  become 
saturated  with  the  water. 

Concrete  is  able  to  retain  water  because  of  the  menisci  formed  by  a wetting  liquid 
with  surface  tension  7 in  pores  of  radius  R.  For  a concrete  of  thickness  L resting  upon  a 
surface  of  water,  the  menisci  within  the  pores  at  the  top  surface  of  the  concrete  specimen 
generate  a pressure  drop  of  27 / R as  one  passes  from  the  air  above  into  the  pore  fluid  just, 
below  the  menisci  surface.  The  pressure  is  sufficient  to  overcome  the  head  pressure  of  pgL , 
thereby  drawing  water  in  from  the  bottom  surface.  This  condition  is  shown  schematically 
in  Fig.  17.  The  pores  in  concrete  are  sufficiently  small  that,  for  concrete  in  contact  with 
water  on  its  bottom  surface  and  in  contact  with  air  at  the  top  surface,  the  concrete 
will  remain  saturated,  and  there  will  be  no  downward  flow  of  pore  fluid.  This  gives  an 
effective  hydraulic  head  pressure  of  zero.  In  fact,  for  an  unsaturated  environment  at  the 
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Figure  18:  Schematic  of  the  hydraulic  pressures  in  a porous  material  with  a thin  layer  of 
water  on  its  top  surface  and  unsaturated  air  exposure  on  its  bottom  surface. 


top  surface,  the  water  will  wick  away  due  to  evaporation,  creating  an  upward  movement 
of  fluid. 


6=3  Break-Through  Pressure 

The  effect  of  the  menisci  can  also  be  applied  to  a concrete  with  water  on  its  top  surface 
and  exposed  to  air  at  the  bottom  surface.  If  the  bottom  surface  of  a concrete  is  exposed 
to  an  unsaturated  environment,  the  bottom  surface  of  the  concrete  will  have  menisci  in 
the  pores,  analogous  to  the  previous  discussion.  Imagine  a very  shallow  pool  of  water 
exists  on  top  of  a concrete  specimen  with  thickness  L.  The  hydraulic  pressure  within  the 
concrete  increases  as  pgz , with  2 represented  the  distance  below  the  surface  of  the  water 
above  the  concrete.  At  the  menisci  on  the  bottom  surface,  there  is  a pressure  increase 
of  approximately  2y / R as  one  passes  from  the  pore  fluid  to  the  air  below  the  specimen. 
This  is  shown  schematically  in  Fig.  18. 

The  minimum  specimen  thickness  L required  for  the  flow  of  water  through  the  spec- 
imen can  be  calculated  from  the  balance  of  these  pressures: 

2y 

P9L  = (45) 

In  SI  units,  p = 1000  kg-m-3,  g = 9.8  m-s-2,  and  7 = 0.0735  N-m-/  (liquid-vapor  surface  at 
15°C).  The  approximate  value  for  the  menisci  radii  corresponding  to  a given  permeability 
can  be  approximated  from  the  Katz-Thompson  relationship: 


k = 


(46) 


The  quantities  k,  F,  and  dc  are  the  specimen  permeability,  formation  factor,  and  critical 
pore  diameter,  respectively.  For  typical  concretes  of  good  quality,  k = 10-17  m2  and 
F = 500.  Using  the  critical  diameter  dc  as  an  estimate  of  the  meniscus  diameter,  this 
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approximation  gives  R — 1.0  nm.  The  corresponding  hydraulic  head  is  15  000  m.  This  is 
also  an  estimate  of  the  hydraulic  head  required  to  overcome  the  effect  of  surface  tension 
or  order  for  How  to  initiate.  And  since  the  hydraulic  pressure  increases  as  pgz  both  inside 
concrete  and  in  free  water,  for  a thin  him  of  water  over  the  concrete,  the  specimen  would 
have  to  be  15  000  m thick  in  order  for  How  to  initiate.  Although  this  calculation  was  a 
rough  approximation,  this  minimum  thickness  is  nearly  four  orders  of  magnitude  greater 
than  one  would  expect  for  a vault  roof.  Therefore,  one  would  not  expect  water  to  How 
through  a vault  roof  that  had  a thin  film  of  water  on  its  top  surface,  and  an  exposure  of 
unsaturated  air  on  the  bottom  surface. 

This  simple  calculation  demonstrates  the  significant  effect  that  an  unsaturated  envi- 
ronment at  the  lower  concrete  surface  can  have  an  advective  transport  through  undam- 
aged concrete.  By  keeping  the  interior  environment  of  the  vault  below  saturation,  the 
advective  flow  through  the  concrete  will  be  virtually  zero. 

6.4  Transport  Through  Cracks 

The  preceding  section  demonstrated  the  significance  of  having  undamaged  concrete  in 
contact  with  a unsaturated  interior  vault.  In  practice,  the  concrete  will  likely  have 
some  distress  that  will  cause  cracking,  and  may  have  a far  different  performance  than 
undamaged  concrete. 

One  can  characterize  the  effects  of  cracks  by  calculating  the  minimum  crack  width 
before  water  will  flow  through  the  concrete.  Consider  a horizontal  concrete  slab  that  is  1 
m thick.  Assume  that  the  slab  is  saturated,  that  there  is  a shallow  film  of  water  on  the 
top  surface,  and  that  the  bottom  surface  is  exposed  to  an  unsaturated  environment  like 
that  shown  in  Fig.  18.  The  pgh  hydraulic  pressure  is  just  able  to  “push”  water  through 
a 15  pm  crack  (assuming  plane,  parallel  walls)  which  would  cause  dripping  from,  or 
pooling  on,  the  bottom  surface  of  the  concrete.  A 15  pin  crack  is  generally  invisible 
to  the  naked  eye,  yet  its  effect  upon  the  advective  transport  through  the  concrete  is 
dramatic.  Clearly,  to  accurately  quantify  advective  transport  through  concrete  one  has 
to  accurately  quantify  the  cracking  within  the  concrete. 

6.4.1  Residence  Time 

Using  this  more  clearly  constructed  description  of  the  effects  of  cracks,  it  should  be  much 
easier  to  characterize  flow  through  the  cracks,  and  to  estimate  a residence  time  for  the 
time  required  for  water  to  traverse  the  concrete  member.  If  there  exist  cracks  either  at 
the  bottom  of  the  slab  or  through  the  entire  slab,  there  should  be  mass  flow,  even  if 
the  atmosphere  of  the  vault  is  unsaturated  with  respect  to  water  vapor.  Knowledge  of 
the  residence  time  is  needed  in  order  to  characterize  the  chemistry  of  the  effluent  into 
the  interior  of  the  vault.  If  the  residence  time  is  sufficiently  short,  the  fluid  will  not 
have  sufficient  time  to  reach  equilibrium  before  passing  through  the  concrete.  This  effect 
directly  effects  any  attempt  to  predict  the  chemistry  of  the  effluent  from  the  vault. 

The  residence  time  for  a crack  can  be  calculated  using  Darcy’s  law.  Within  a vertical 
crack,  there  is  a body  force  of  pg  downward  on  the  water.  This  corresponds  to  a constant 


43 


pressure  gradient:  VP  = —pg.  Darcy’s  law  gives  the  bulk  fluid  velocity  v6: 

v„  = — — VP 
P 


(47) 


The  averaged  pore  velocity  vp,  is  greater  than  the  bulk  velocity  by  a factor  of  the  porosity: 


v 


p 


V- 

0 


(48) 


Therefore,  the  residence  time  tr  for  a fluid  to  pass  through  a sample  of  thickness  L is 
tr  — L/vp.  For  a quality  concrete  ( k — 10-18  m2)  1 m in  thickness,  with  1 m of  head, 
the  Darcy  velocity  is  less  than  1 p,m/d.  The  corresponding  residence  time  is  2 700  years. 

For  cracked  concrete,  the  Darcy  velocity  is  much  greater.  As  an  example,  consider 
a 1 m thick  slab  with  a 100  pm  crack  extending  through  the  entire  slab.  Since  the 
permeability  of  plane,  parallel  surfaces  separated  by  a distance  w is  w2)  12  [47],  the 
permeability  of  the  100  pm  crack  is  8.3xlO-10  nr.  Assuming  the  pore  fluid  density  and 
viscosity  are  approximately  equal  to  those  of  water,  for  an  equivalent  water  head  of  1 m, 
the  bulk  velocity  vb  through  the  crack  is  8.3  mm/s.  For  the  100  pm  crack,  vp  - vb,  and 
the  residence  time  tT  is  approximately  120  s.  Therefore,  the  issue  of  pore  fluid  residence 
time  in  cracked  concrete  is  important. 


6.4.2  Evaporative  Mass  Transfer 


In  cases  where  a concrete  is  kept  saturated  and  the  there  are  cracks  at  on  the  interior 
surface,  water  will  probably  advance  to  the  surface  of  the  concrete  where  these  cracks 
exist.  Given  the  size  of  typical  cracks,  the  equilibrium  relative  humidity  of  an  atmo- 
sphere exposed  will  be  nearly  100  % RH.  Therefore,  a vault  interior  environment  below 
saturation  will  cause  evaporation  of  the  concrete  pore  fluid  that  is  in  the  cracks. 

Assuming  that  the  concrete  permeability  of  the  uncracked  portion  is  sufficiently  great 
so  as  to  supply  these  cracks  with  additional  pore  solution,  there  should  be  a thin  film  of 
water  on  the  interior  surface  of  the  concrete.  The  rate  of  evaporation  from  the  bottom 
surface  can  be  approximated  from  formulas  for  evaporation  from  a free  surface  of  water. 
The  saturation  pressure  Ps  (Pa)  of  water  can  be  expressed  as  a function  of  its  temperature 
T (in  Celsius)  [46]: 


Ps{T)  — 610  exp 


• 17.3  T • 
.237.3  + T. 


(49) 


Let  Ps(Tc)  and  Ps(Ta)  represent  the  saturation  pressure  of  water  at  the  temperature 
of  the  concrete  (Tc)  and  of  the  air  (Ta),  respectively.  Exposed  to  an  atmosphere  with 
relative  humidity  £ and  with  wind  speed  (forced  convection,  if  applicable)  vw,  (m  s-1) 
across  the  concrete  surface,  the  water  evaporation  rate  E (kg-m~2-s_1)  is  [46] 


E = (8.69  x 10~5)  (PJ  - £P£)(0.253  + 0.017  vw)  (50) 

For  an  underground  concrete  structure,  the  temperature  of  the  concrete  would  be  equal 
to  the  vault  temperature.  Therefore,  the  evaporation  rate  equation  for  a concrete  vault 
at  temperature  T simplifies  slightly: 

E(T)  = (8.69  x 10~5)  PS{T)  (1  - f)  (0.253  + 0.017uw)  (51) 
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This  is  the  most  general  result.  If  the  interior  atmosphere  of  the  vault  is  maintained 
constant  through  convection,  the  wind  speed  will  be  non-zero.  If  the  vault  interior  is 
sealed,  the  evaporation  rate  will  reach  zero  over  time. 

Once  the  evaporation  rate  is  determined,  it  can  be  compared  with  the  advective  flux. 
If  the  evaporation  rate  is  less  than  the  advective  flux,  there  will  exist  accumulations  of 
pore  solution  on  the  bottom  surface  of  the  roof,  and  this  accumulation  will  drip  onto  the 
contents  of  the  vault. 
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7 Predicting  Drying  Shrinkage  and  Flexural  Cracks 

A recent  review  of  the  4SIGHT  computer  model  [47]  revealed  that  the  most  important 
factor  affecting  ionic  transport  through  concrete  is  the  presence  of  cracks.  The  presence 
of  cracks  not  only  controls  the  quantity  of  ions  transported,  but  can  also  control  whether 
there  will  be  convective  transport  in  the  absence  of  externally  applied  hydraulic  pressure. 
That  recent  review  is  reproduced  here  nearly  in  its  entirety  for  completeness. 

Cracking  can  originate  from  within  the  concrete  due  to  a number  of  mechanisms. 
During  the  placement  of  the  concrete,  if  the  evaporation  rate  is  great  enough,  the  concrete 
surface  can  develop  tensile  stresses  sufficient  to  crack  the  plastic  concrete.  These  plastic 
shrinkage  cracks  typically  extend  through  the  section  of  the  concrete.  If  during  placement 
the  concrete  member  experiences  settlement,  and  the  concrete  has  insufficient  strength 
to  withstand  the  resulting  stresses,  the  concrete  will  undergo  settlement  cracking.  If  the 
concrete  member  is  a slab  that  spans  a distance,  the  concrete  surface  in  tension  (typically 
the  bottom  surface)  will  possibly  develop  flexural  cracks  after  removal  of  the  concrete 
form  work.  During  the  cement  paste  hydration,  a large  portion  of  the  initial  mix  water 
will  be  consumed,  with  the  generation  of  heat.  Differences  between  concrete  temperature 
and  ambient  temperature  can  lead  to  thermal  cracking , depending  on  the  mechanical  state 
of  the  concrete.  The  continued  removal  of  water  by  the  hydration  process  will  generate 
a chemical  shrinkage  stress  that  can  initiate  autogeneous  shrinkage  cracks.  Subsequent 
drying  due  to  ambient  conditions  will  also  generate  a shrinkage  strain  that  can  initiate 
drying  shrinkage  cracks  when  the  concrete  is  restrained. 

Some  cracks  are  to  be  expected,  while  others  are  the  result  of  improper  construction 
practices.  Plastic  shrinkage  cracking,  thermal  cracking,  and  settlement  cracking  can  be 
either  avoided  or  minimized  by  following  industry  guidelines  (ACI  Manual  of  Concrete 
Practice)  for  the  placement  and  curing  of  concrete.  Both  flexural  cracks  and  shrinkage 
cracks  are  to  be  expected  (under  the  assumptions  given  above),  and  the  resulting  crack 
widths  and  spacings  can  be  estimated.  Therefore,  the  incorporation  of  crack  modeling 
by  4SIGHT  will  be  limited  to  the  prediction  of  flexural  and  shrinkage  cracking  and  their 
effect  on  transport. 

The  models  for  the  crack  width  and  spacing  described  here  will  be  implemented 
into  a structural  analysis  of  a vault  that  is  meant  to  represent  a typical  design.  Under 
ideal  conditions,  a structural  engineer  would  perform  an  analysis  of  the  specific  vault  in 
question,  and  then  input  the  calculated  parameters  into  4SIGHT.  The  structural  analysis 
presented  is  not  meant  to  represent  an  exhaustive  analysis.  It  is  the  default  analysis, 
however,  that  will  be  performed  in  the  absence  of  any  analysis  by  the  engineer  who  is 
using  4SIGHT. 

To  proceed,  one  must  assume  an  approximate  geometry  for  a representative  vault 
structure.  The  principles  used  here  for  estimating  the  composite  permeability  due  to 
cracking  are  general  enough  to  apply  to  many  different  vault  designs.  A general  geometry 
and  construction  plan  for  the  vault  will  be  assumed  here  to  illustrate  the  analysis. 

To  partially  simplify  the  structural  analysis,  the  vault  will  be  approximated  by  a 
container  with  horizontal  dimensions  L i and  L2  such  that  their  ratio  is  greater  than 
two:  Lx/Lo  > 2.  The  corresponding  edge  loads  / 1 and  /2  have  a quartic  dependence  on 
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length  [48]: 


.f±  = IA 

h L\ 


(52) 


Since  the  load  carried  in  the  direction  of  the  shorter  span  is  16  times  greater,  the  vault 
analysis  will  be  performed  using  one-way  slab  behavior. 

The  spacing  and  width  of  both  shrinkage  and  flexural  cracks  depends  strongly  on  the 
amount  and  the  location  of  the  steel  reinforcement.  In  addition,  the  positive  moment 
due  to  the  load  will  depend  on  how  the  roof  is  attached  to  the  vault  walls.  The  roof 
could  either  rest  on  the  walls  (simply  supported)  or  it  could  be  rigidly  attached  to  the 
walls  (continuous  reinforcement).  For  a given  geometry  and  load,  the  simply  supported 
roof  will  subjected  to  approximately  twice  the  positive  moment  of  a fixed  end  roof. 

Although  an  analysis  that  assumes  a simply  supported  roof  would  result  in  a greater 
estimated  positive  moment,  it  is  quite  unlikely  that  simply  supported  construction  would 
occur  in  practice.  By  definition,  a simply  supported  roof  must  be  able  to  rotate  and 
translate  freely  over  the  vault  walls.  While  this  type  of  construction  can  eliminate  drying 
shrinkage  cracks,  the  joint  provides  a likely  path  for  the  ingress  of  fluids  and  ionic  species. 

In  addition  to  the  analysis  of  the  flexural  and  shrinkage  crack  formation  at  the  center 
of  the  roof,  there  will  exist  negative  moments  at  the  connection  between  the  roof  and 
the  walls.  Although  an  analysis  of  negative  moment  cracks  is  not  performed  here,  the 
method  for  incorporating  these  cracked  sections  into  the  overall  composite  properties 
follows  from  the  positive  moment  analysis.  In  addition,  no  analysis  of  the  shear  stress  at 
the  joint  between  the  roof  and  the  wall  is  presented. 

The  scenario  adopted  here  is  that  the  vault  will  first  be  filled  with  waste  materials, 
compacted,  and  the  roof  cast  in  place  on  top  of  the  vault  contents.  The  expected  thickness 
of  the  vault  roof  is  assumed  to  be  1 m,  and  the  contents  of  the  vault  will  probably  settle 
under  the  dead  load  of  the  roof  and  the  soil  above  the  vault.  The  roof  will  then  flex, 
possibly  forming  flexural  cracks.  In  fact,  a small  amount  of  flexure  will  be  useful  in  that 
it  will  help  close  shrinkage  cracks  at  the  top  surface. 


7.1  Cracked  Reinforced  Concrete 

Given  the  construction  plan  outlined  above,  the  concrete  roof  will  likely  experience  shrink- 
age stresses  during  the  first  months,  and  later  experience  flexural  stresses  as  the  vault 
contents  settle  over  time.  While  the  flexural  stresses  will  widen  the  shrinkage  cracks  on 
the  bottom  of  the  roof,  they  will  contract  the  shrinkage  cracks  on  the  top  surface  of  the 
vault.  The  reinforced  concrete  containing  both  shrinkage  and  flexural  cracks  will  have  a 
permeability  that  is  greater  than  the  uncracked  concrete. 

7.1.1  Crack  Model 

For  simplicity,  the  cracks  are  approximated  by  plane  parallel  walls  separated  by  the 
observed  crack  width.  The  open  cracks  on  the  tension  face  of  the  roof  propagate  upwards 
to  the  neutral  axis.  While  the  depth  of  crack  penetration  is  accurate,  the  approximation 
of  the  cracks  by  plane  parallel  walls  is  conservative.  This  assumption  does  not  account 


for  either  the  “V”  shape  of  the  crack  or  the  crack  wall  surface  roughness,  both  of  which 
would  act  to  reduce  the  crack  permeability  with  respect  to  one  with  plane  parallel  sides. 

The  permeability  k'  of  a crack  with  plane  parallel  sides  is  solely  a function  of  the 
crack  width  w [47]: 
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7.1.2  Composite  Permeability 
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Figure  19:  Composite  models  for  the  bulk  permeability  of  the  roof  incorporating  the 
uncracked  permeability  kc,  crack  permeability  k',  and  having  a number  m cracks  over 
the  slab  length  L:  (a)  slabs  with  m cracks,  each  with  width  w\  (b)  slabs  with  m\  and  m2 
cracks  of  width  u)\  and  W2,  respectively. 


Composite  models  for  the  roof  bulk  permeability  are  shown  in  Fig.  19.  The  cracks 
propagate  from  the  bottom  of  the  slab  up  to  the  neutral  axis.  The  thickness  of  the  slab 
is  h,  and  the  depth  to  the  neutral  axis  is  expressed  as  the  product  ah.  The  permeability 
of  the  uncracked  concrete  is  kc.  The  composite  has  an  additional  component  with  per- 
meability k'  due  to  m cracks,  each  with  width  w,  over  the  entire  span  L.  For  Fig.  19(a), 
the  permeability  kb  of  the  two  bottom  sections  kc  and  k'  in  parallel  is  a weighted  average: 
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The  composite  permeability  k*  of  kb  and  kc  (top  section  in  Fig.  19(a)  ) in  series  is  treated 
like  the  corresponding  problem  of  conductors  in  series: 
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Given  that  a crack  with  a width  of  one  micrometer  has  a permeability  of  8x10  1!  nr, 
and  that  concrete  permeability  is  on  the  order  of  10  18  nr,  the  second  term  in  Eqn.  56 
can  be  neglected: 

k*(ki>  » kc)  = — (56) 

a 

Therefore,  for  a given  concrete  mixture,  the  control  of  roof  permeability  is  the  control 
of  the  neutral  axis  depth,  given  that  there  is  sufficient  flexural  strain  to  close  shrinkage 
cracks  at  the  top  surface  of  the  roof. 

The  relationship  between  composite  permeability  and  the  neutral  axis  depth  should 
not  be  used  as  a means  to  reduce  vault  permeability.  Within  practical  limits,  the  depth 
of  the  neutral  axis  can  only  affect  the  composite  permeability  by  a factor  of  two  or 
three.  However,  variations  in  concrete  mixture  proportions  and  cementitious  materials 
can  change  the  concrete  permeability  by  a factor  of  10  to  100.  Therefore,  efforts  to  reduce 
overall  permeability  should  concentrate  on  controlling  the  concrete  mixture. 

7.1.3  Composite  Permeability  Model 

The  way  in  which  the  effects  of  shrinkage  and  flexural  cracks  are  combined  depends  upon 
the  formation  and  the  geometry  of  the  cracks.  The  shrinkage  cracks  will  likely  form 
first,  and  there  will  be  a number  ms  of  these  cracks,  each  with  width  ius,  separated  by 
a distance  Ls.  Similarly  for  the  flexural  cracks,  formed  independently  of  the  shrinkage 
cracks,  there  will  be  m j flexural  cracks,  each  with  width  wj , separated  by  a distance 
Lf.  The  objective  is  to  consider  the  formation  of  the  cracks  independently,  and  then 
determine  the  properties  of  the  cracked  material. 

The  result  ing  width  and  separation  of  cracks  depends  on  the  cracking  scenario.  Here, 
the  assumption  is  that  the  shrinkage  cracks  form  first  and  the  flexural  cracks  form  after 
the  vault  contents  settle.  For  the  case  when  the  crack  spacings  are  approximately  equal, 
so  that  Ls  ~ Lf , the  later  flexural  cracks  should  simply  widen  the  existing  shrinkage 
cracks.  If  Ls  <C  T/,  the  flexure  should  once  again  widen  the  existing  shrinkage  cracks. 
For  the  third  case  when  Ls  Lf , newr  flexural  cracks  should  form  in  between  the 
existing  shrinkage  cracks,  each  with  width  Wf  and  separated  by  Lf.  If  one  assumes  that 
the  existing  shrinkage  cracks  contribute  to  the  flexural  cracking,  the  shrinkage  cracks 
should  also  widen  by  an  amount  equal  to  Wf. 

The  particular  crack  model  to  choose  from  in  Fig.  19  depends  upon  the  ratio  LJLf. 
For  the  first  two  scenarios  presented  ( Ls/Lf  ~ \,LS/Lf  <C  1),  the  flexural  stress  widens 
the  existing  shrinkage  cracks.  Each  crack  will  then  have  the  same  width  and  spacing, 
corresponding  to  the  model  shown  in  Fig.  19(a).  For  the  third  case  ( Ls/Lf  ^>1),  there 
will  exist  two  crack  widths  and  spacings.  The  existing  shrinkage  cracks  will  widen  by  an 
amount  Wf,  and  the  nif  — ms  flexural  cracks  will  each  have  a width  Wf,  corresponding  to 
the  model  shown  in  Fig.  19(b).  The  analysis  of  this  ternary  composite  bulk  permeability 
k*  follows  from  the  binary  example  given  previously.  The  parameters  for  each  model  are 
shown  in  Table  9. 
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Table  9:  Parameters  for  the  composite  permeability  calculation  for  m cracks  of  width  w. 
The  crack  model  refers  to  that  shown  in  either  Fig.  19(a)  or  Fig.  19(b). 
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7.2  Flexural  Crack  Model 

The  model  for  estimating  flexural  crack  widths  and  spacings  is  based  upon  the  recom- 
mendations given  in  ACI  224.2R-90  [3].  The  recommendation  for  estimating  flexural 
crack  width  wj  is  the  Gergely-Lutz  equation: 

wf  = 2.2  0er(deAn)1/3  (57) 


The  definitions  of  the  parameters  are  as  follows: 


d.\  : distance  from  the  neutral  axis  to  the  tension  face 

g?2  •’  distance  from  the  neutral  axis  to  the  centroid  of  the  tension  steel  reinforcement 
ft  : the  ratio  d\/d2 

er  : strain  in  the  tension  steel  reinforcement 

dc  : minimum  concrete  cover  thickness,  measured  to  center  of  tension  reinforcement 
An  : the  largest  area  of  concrete,  perpendicular  to  the  reinforcement,  sym- 
metric about  the  centroid  of  the  tension  reinforcement,  divided  by  the 
number  of  bars 


The  estimate  of  the  spacing  between  flexural  cracks  is  based  upon  the  CEB/FIP 
model.  The  maximum  distance  ls  over  which  slip  can  occur  can  be  approximated  from 
the  following  equation: 


Is 


dr 

3.6  ps 


(58) 


The  quantity  dr  is  the  steel  reinforcement  diameter,  and  ps  is  the  area  of  steel  divided 
by  the  effective  area  of  concrete  in  tension.  The  average  spacing  between  flexural  cracks 
Lf  is  taken  to  be  2/3  this  value: 


Lf  = 


dT 

5-4  ps 


(59) 


7.3  Shrinkage  Crack  Model 

The  drying  shrinkage  model  is  based  upon  the  calculations  of  Base  and  Murray  [49]  for 
calculating  the  width  ws  and  spacing  Ls  of  shrinkage  cracks.  The  model  assumes  the 
member  is  restrained  at  its  ends  and  requires  the  following  design  parameters: 
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L : length  of  slab 

Es  : Young’s  modulus  of  the  steel  reinforcement 
Ec  : Young’s  modulus  of  the  concrete 
n : modulus  ratio  = Es/Ec 
As  : total  area  of  steel  in  cross  section 
Ac  : total  area  of  concrete  in  cross  section 
b : slab  width 
d : depth  of  reinforcement 
p : reinforcement  ratio  = As/bd 
dT  : steel  reinforcement  bar  diameter 
esh  : shrinkage  strain 

et  : concrete  tensile  strain  at  crack  initiation 


The  “no-bond”  or  slip  length  a of  reinforcement  near  a shrinkage  crack  can  be  ap- 
proximated from  the  reinforcement  diameter  dr  and  the  reinforcement  ratio  p: 


dT 

a = 0.08  — 
V 


(60) 


The  number  of  shrinkage  cracks  rns  that  form  along  a restrained  slab  of  length  L can  be 
estimated  from  the  following  equation: 


ms  = 1 + 


Lap 

~2a 


fs/l 

- 3 et 
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The  spacing  between  shrinkage  cracks  Ls  can  be  estimated  from  the  number  of  shrinkage 
cracks  ms  over  the  slab  span  L: 


Ls 


L 

ms  + 1 


(62) 


The  estimate  for  the  number  of  shrinkage  cracks  ms  can  be  used  to  estimate  the  steel 
stress  fs: 


L — 2 ma  \ 
npL  + 2 ma  ) 


(63) 


The  average  shrinkage  crack  width  ws  can  be  estimated  from  the  following  equation: 


ws  = 2 a 


fs  ^sh  \ 

E,  ~3fJ 


(64) 


7.4  Design  Example 

The  following  is  an  example  of  how  one  can  perform  a calculation  to  estimate  the  spacing, 
width,  and  depth  of  shrinkage  and  flexural  cracks.  The  calculation  begins  with  the 
concrete  mixture  design,  estimates  the  extent  of  cracking,  and  then  estimates  of  the 
composite  permeability. 
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7.4.1  Mixture  Design 

The  mixture  design  chosen  for  this  example  represents  a balance  between  quality  of 
concrete  and  ease  of  formulation  and  placement.  The  design  compressive  strength  /'  is 
35  MPa,  which  is  within  the  normal  strength  guidelines  of  ACI  211. 

The  constituent  material  properties  are  as  follows: 


Cement  Type 
Cement  Density 
CA  Density 
FA  Density 
Max.  Agg.  Diameter 


Type  I 
3.2  g/cm3 
2.8  g/cm3 
2.6  g/cm3 
25  mm 


For  a reinforced  foundation  with  a design  strength  of  35  MPa,  the  ACI  211  procedure 
results  in  the  following  trial  mixture  proportion: 


Constituent  Mass  (kg/m3) 


Water 

190 

Cement 

399 

Coarse  Agg. 

1104 

Fine  Agg. 

757 

A high  range  water  reducer  is  not  specified,  simplifying  the  formulation. 

Although  this  mixture  is  only  a trial  proportion,  it  will  be  used  as  the  final  mixture 
proportion  for  the  purposes  of  demonstration. 


7.4.2  Concrete  Properties 

Based  on  the  ACI  211  procedure,  the  concrete  is  assumed  to  have  the  following  properties: 

w/c  0.43 

Slump  75  mm 

Air  Content  1.5  % 

Concrete  Density  2450  kg/m3 

In  practice,  these  properties  would  be  determined  from  tests  performed  on  trial  batches 
and  not  based  on  the  ACI  211  procedure.  As  a demonstration,  however,  the  ACI  values 
will  be  used  here. 

The  equation  for  estimating  the  elastic  modulus  of  concrete  is  described  in  Nilson  [48]: 


Ec  = ^3320  [MPa1/2]///  + 6890  [MPa] 


Pc 


1.5 


2330  [kg/nr 


(65) 


This  equations  can  be  applied  to  normal  density  concretes  over  a wide  range  of  design 
strengths:  21  MPa  < f'c  < 83  MPa.  For  the  mixture  given  above,  the  estimated  elastic 
modulus  is  28.6  GPa. 
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A concrete  member  in  tension  will  experience  a tensile  stress  ft.  In  flexure,  the 
member  will  remain  uncracked  until  the  tensile  stresses  exceed  the  modulus  of  rupture. 
The  estimation  for  the  modulus  of  rupture  fr  is  taken  from  Nilson  [48]: 

0.67  [MPa1'2]//2  <fr<  1.0  [MPa1/2]//2 

For  the  mixture  used  in  this  example,  it  is  expected  that  the  modulus  of  rupture  fr  can 
be  characterized  by  the  following  inequality: 

4.0  MPa  < fr  < 5.9  MPa 


Figure  20:  Quantities  for  estimating  shrinkage  strain:  (a)  ultimate  drying  shrinkage 
strain  under  standard  conditions;  (b)  factor  ycp  for  duration  of  moist  curing;  (c)  factor 
70  for  ambient  relative  humidity;  and  (d)  factor  7„s  for  the  volume  to  surface  ratio. 


The  ACI  guideline  for  estimating  the  ultimate  drying  shrinkage  strain  eshM  is  shown 
in  Fig.  20(a)  as  a function  of  the  mixture  water  content.  For  this  mixture,  the  expected 
ultimate  drying  shrinkage  is  approximately  600  gm/m. 
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The  expected  shrinkage  strain  esh  will  be  less  than  the  ultimate  value  due  to  mitigating 
factors.  ACI  209R-82  gives  the  following  mathematical  expression  for  the  combined 
effects  of  these  factors  on  the  expected  ultimate  shrinkage  strain: 

C sh  £sh,u  J[  T (66) 


The  quantity  eshiU  is  the  ultimate  shrinkage  strain  under  standard  conditions:  curing 
duration  of  7 d,  exposure  to  40  % relative  humidity,  and  a volume-to-surface  ratio  of 
0.038  m. 

There  are  three  factors  y,  considered  here:  duration  of  curing  ycp;  ambient  relative 
humidity  y^;  and  volume-to-surface  ratio  jvs.  There  are  other  factors  outlined  in  ACI 
209,  but  the  factors  relating  to  concrete  mixture  design  do  not  apply  here  since  we  are 
using  the  ultimate  drying  shrinkage  prediction  shown  in  Fig.  20(a). 

The  factor  ycp  for  the  duration  of  moist  curing  is  shown  in  Fig.  20(b)  as  a function  of 
curing  time  t.  The  filled  circles  are  the  recommendations  from  ACI  209,  and  the  curve 
is  an  approximation: 

7cp  = 1 - 0.103  In  ^ (67) 

The  value  of  r is  7 d,  corresponding  to  the  standard  curing  condition  for  the  duration  of 
moist  curing. 

The  factor  y^  for  the  ambient  environment  is  shown  in  Fig.  20(c)  as  a function  relative 
humidity  0.  Over  the  range  of  relative  humidity  shown  in  Fig.  20(c),  the  factor  has  two 
linear  regions: 

f 1.4  - 1.00  0.40  < 0 < 0.80 

= { (68) 

[ 3.0  — 3.0  0 0.80  <4><  1.00 

It  is  difficult  to  estimate  the  ambient  relative  humidity  for  a buried  vault.  However,  most 
certainly  the  relative  humidity  will  be  greater  than  40  % and  less  than  100  %.  As  a 
conservative  estimate,  the  factor  y^  will  be  assigned  a value  of  one  here.  However,  in 
practice,  an  assessment  should  use  a value  commensurate  with  the  specific  soil  conditions 
being  considered. 

The  factor  yvs  for  the  geometry  of  the  structural  element  is  shown  in  Fig.  20(d)  as  a 
function  of  the  volume-to-surface  ratio  v.  ACI  209  expresses  this  factor  in  analytic  form: 


7vs  = F2  exp  {-vjv0) 


(69) 


The  quantity  v0  equals  0.212  m,  and  simply  renormalizes  v.  For  the  present  design,  the 
length  and  width  of  the  slab  are  L and  2 L,  respectively.  For  a slab  thickness  h,  the 
volume-to-surface  ratio  v is  the  following: 


h 

2 + 3f 


(70) 


For  a slab  thickness  h equal  to  1 m,  and  a slab  length  L equal  to  10  m,  the  volume-to- 
surface  ratio  v equals  0.435  m,  corresponding  to  a shrinkage  factor  yvs  of  0.15. 
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The  expected  ultimate  shrinkage  strain  esb  is  the  product  of  the  ultimate  shrinkage 
strain  under  standard  conditions  esh,u  and  the  factors  7Z: 


f sh  ^sh,u  'icp  Ivs 

— 90  fim / in 


(71) 


This  is  significantly  less  than  the  GOO  //m/m  that  was  based  on  standard  conditions. 

7.4.3  Steel  Properties 

The  steel  considered  for  this  example  is  ASTM  A 615  Grade  60  steel;  a steel  with  a 
60  000  psi  yield  strength.  The  corresponding  properties  in  SI  units  are  as  follows: 


Es  = 200  GPa 
fy  = 414  MPa 


(72) 


Therefore,  the  modulus  ratio  n,  defined  as  Es/Ec , is  approximately  7. 

7.4.4  Reinforcement  Ratio  Requirements 

One  of  the  most  important  parameters  in  estimating  the  bulk  permeability  of  the  re- 
inforced concrete  element  is  the  depth  of  the  neutral  axis.  This  quantity  is  estimated 
from  a structural  analysis  of  the  complete  concrete  structure.  Ideally,  an  engineer  using 
4SIGHT  will  have  access  to  a structural  analysis  that  includes  an  estimate  of  the  depth 
of  the  neutral  axis.  In  the  absence  of  this  information,  4SIGHT  must  make  a reasonable 
approximation. 

The  following  section  is  not  meant  as  a primer  on  structural  analysis.  Rather,  it 
outlines  the  procedure  4SIGHT  will  use  to  estimate  the  neutral  axis  in  the  absence  of 
such  data.  If  the  engineer  using  4SIGHT  determines  that  the  analysis  invoked  by  4SIGHT 
is  insufficient,  the  engineer  must  then  perform  a more  thorough  analysis  independently, 
and  report  the  neutral  axis  as  input  to  the  4SIGHT  computer  program. 

There  is  a limit  on  both  the  minimum  and  maximum  amount  of  steel  reinforcement 
allowed  in  the  design  of  a flexural  member.  The  minimum  limit  is  based  on  the  ability 
of  the  element  to  carry  the  intended  load.  The  maximum  limit  is  determined  by  the 
ultimate  failure  mode  within  the  reinforced  concrete  element. 

The  two  regions  of  stress  within  the  slab  are  the  tension  zone  in  the  bottom  and  the 
compression  zone  in  the  top.  At  the  ultimate  load,  either  the  steel  will  fail  in  tension,  or 
the  concrete  will  fail  in  compression.  Steel  failure  in  tension  is  a gradual  process,  giving 
“warning”  of  impending  failure.  Compressive  failure  of  concrete,  by  contrast,  is  a sudden 
mode  of  failure.  Therefore,  there  is  a maximum  limit  on  the  amount  of  steel  reinforcement 
so  as  to  ensure  that  the  steel  will  yield  before  the  concrete  fails  in  compression. 

The  amount  of  steel  at  which  the  two  modes  of  failure  coincide  at  ultimate  load  is 
the  balanced  reinforcement  ratio  pb: 


(73) 
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For  a 35  MPa  concrete  mixture,  the  recommended  value  of  is  approximately  0.775. 

The  maximum  allowable  reinforcement  ratio  pmax  is  3/4  of  the  balance  reinforcement 
ratio: 


Pmax  0.75  pi) 
= 0.025 


(74) 


The  minimum  reinforcement  ratio  pmm  can  be  calculated  based  on  the  assumption 
that  a steel  stress  fs  that  is  2/3  of  the  yield  stress  is  an  acceptable  limit  of  steel  stress 
immediately  after  cracking.  The  corresponding  reinforcement  is  the  following: 

0.220  [MPa'/2] 

Pmin  r \ h) 

Jy 

= 0.0031 


This  is  in  agreement  with  the  ACI  318  recommendation: 


Pmin 


1.38  [MPa] 

~/T 

= 0.0033 


(76) 


7.4.5  Roof  Geometry 


(a)  (b)  (c) 

Figure  21:  Roof  geometry  used  for  analysis:  (a)  cross  section  showing  location  of  steel 
reinforcement  (filled  circles);  (b)  transformed  section  used  for  uncracked  analysis;  (c) 
transformed  section  used  for  cracked  analysis. 

The  dimensions  of  the  roof  slab  are  defined  in  Fig.  21  and  are  given  the  following 
values: 
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h : 1.0  m 
d : 0.9  m 
b : 1.0  m 
dr  : 35.8  mm 

The  slab  width  b is  set  to  1.0  m to  simplify  the  analysis  by  normalizing  values  on  a “per 
meter”  basis.  The  reinforcement  diameter  dr  of  35.8  mm  corresponds  to  a No.  36  metric 
bar  size  or  a No.  11  bar  in  customary  inch-pound  units. 

7.4.6  Positive  Moment 

The  positive  moment  at  the  center  of  the  roof  slab  is  determined  from  the  dead  load 
due  to  the  concrete  mass  and  the  mass  of  soil  above  the  vault.  The  thickness  h of  the 
concrete  roof  is  1 m,  and  the  assumed  thickness  of  soil  hs  over  the  roof  is  10  m.  The 
density  of  the  concrete  £c  is  2450  kg/nr?,  and  the  density  of  the  soil  £s  is  assumed  to  be 
1500  kg/m,!  (this  soil  density  value  is  used  for  demonstration  purposes  and  is  not  meant 
to  be  an  accurate  estimate  of  any  particular  type  of  soil).  The  nominal  uniform  load 
(pressure)  wn  on  the  roof  depends  upon  the  thickness  of  each  element: 

wn  = tcgh-\-£sghs  (77) 

= 171, 000  N/m2 

The  quantity  g is  the  gravitational  acceleration  (9.8  m/s2).  The  design  load  wu\ 

wu  = 1.4  wn  (78) 

= 239, 000  N/m2 

The  quantity  1.4  is  the  load  factor  for  dead  loads. 

The  positive  moment  M+  at  mid-span  due  to  the  ultimate  load  wu  is  the  integral  of 
the  load  over  half  the  span: 


M+  = 


The  quantity  6 is  a factor  that  depends  on  the  end  conditions  of  the  slab.  For  a simply 
supported  one-way  slab,  6 = 1.  Here,  it  will  be  assumed  that  the  slab  edge  condition 
corresponds  to  the  end  slab  of  a multi-span  roof;  one  end  is  restrained  and  the  other  is 
simply  supported.  Under  these  conditions,  6 — 4/7: 

M+  = 1,  707,  000  N-m 


rb  rL/2 

0 / lwu  dl  db 

Jo  Jo 


0b 


1 0 ./0 
T 9 

wq,L 


(79) 


7.4.7  Uncracked  Analysis 

The  starting  point  for  the  structural  analysis  is  an  analysis  of  the  roof  assuming  that  there 
are  neither  flexural  nor  shrinkage  cracks.  The  maximum  flexural  stress  from  the  design 
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loads  can  then  be  compared  to  the  modulus  of  rupture.  If  the  flexural  stress  exceeds  the 
modulus  of  rupture,  the  analysis  must  be  repeated  assuming  flexural  cracking. 

As  the  slab  flexes,  the  top  portion  of  the  roof  will  be  in  compression,  and  the  bottom 
portion  in  tension.  There  will  exist  a 2-D  surface  through  the  midsection  of  the  slab 
that  defines  where  the  concrete  is  neither  in  tension  nor  in  compression.  This  surface  is 
referred  to  as  the  neutral  axis. 

For  a member  in  pure  bending  (no  shearing  force)  the  strains  are  linearly  proportional 
to  the  distance  from  the  neutral  axis.  Since  the  coefficient  of  proportionality  is  same  for 
both  the  tensile  and  the  compressive  strains  (within  the  elastic  limit),  the  location  of  the 
neutral  axis  can  be  determined  by  balancing  the  area  moments  of  the  slab  cross  section. 

To  simplify  the  calculation  of  the  tensile  and  compressive  area  moments,  the  reinforced 
concrete  slab  cross  section  is  first  transformed  into  an  equivalent  concrete  cross  section. 
Specifically,  the  steel  is  replaced  mathematically  by  an  area  of  concrete  equal  to  (n—  1)AS, 
and  located  at  the  depth  of  the  reinforcement.  This  is  shown  schematically  in  Fig.  ‘21(b) 
for  this  analysis. 

The  equation  for  calculating  the  area  moment  is  the  integral  of  the  moment  arm  r 
(absolute  distance  from  the  neutral  axis)  over  the  slab  cross  section  (da  = dr  db): 


M 


r da 


(80) 


The  moment  of  the  area  above  the  neutral  axis,  which  is  located  at  depth  c from  the  top 
of  the  slab,  is  set  equal  to  the  moment  of  the  area  below  the  neutral  axis: 

be2  b(h  — c )2 

— = ■ 2 + (d-  c)  (n  - 1)  A5  (81) 

The  quadratic  terms  in  c cancel,  yielding  an  algebraic  solution: 

bh 2 + 2(n  — 1 )Asd 


c — 


2bh  + 2 (n  — 1 )AS 

The  steel  area  can  be  expressed  as  a function  of  the  steel  reinforcement  ratio  p: 

A, 


P = 


bd 


This  can  be  used  to  simplify  the  equation  for  the  neutral  axis  depth  c: 


d 


c = 


(^)2  + ‘2(n-l)p 
h/d  + (n  — 1 )p 


(82) 


(83) 


(84) 


The  neutral  axis  depth  c,  as  a function  of  the  reinforcement  ratio  p,  is  shown  in  Fig.  22. 
For  the  uncracked  analysis,  the  depth  of  the  neutral  axis  is  insensitive  to  the  reinforcement 
ratio. 

The  stresses  within  the  transformed  concrete  section  will  depend  upon  the  moment 
of  inertia  / of  the  transformed  slab.  The  moment  of  inertia  about  the  neutral  axis  is  the 
integral  of  the  square  of  the  moment,  arm  over  the  slab  cross  section: 

I — J r2  da  (85) 
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Figure  22:  Depth  of  neutral  axis  c as  a function  of  the  reinforcement  ratio  p for  the 
uncracked  analysis. 


For  the  transformed  section,  the  moment  of  inertia  is 


b(h  - y )3 
3 


+ (d-y)2  {n 


1 )AS 


(86) 


For  linear  elastic  behavior,  the  compressive  stress  fc  and  the  tensile  stress  ft  at  the 
top  and  bottom  surfaces  are  proportional  to  the  distances  from  the  neutral  axis: 


fc  = M+  J f,  = M+  b-pt  (87) 

These  quantities  are  shown  in  Fig.  23  as  a function  of  the  reinforcement  ratio  p.  The 
compressive  stress  fc  is  approximately  one-quarter  the  design  strength  f'c  for  all  values 
of  the  reinforcement  ratio.  By  contrast,  the  tensile  stress  ft  is  considerably  greater  than 
the  estimated  modulus  of  rupture  range  of  4 MPa  to  6 MPa.  Therefore,  one  expects  that 
the  roof  slab  would  develop  cracks  in  the  middle  portion  of  the  slab. 


7.4.8  Cracked  Section  Analysis 

Based  upon  the  previous  analysis,  it  is  likely  that  the  concrete  roof  slab  will  crack  due  to 
flexural  loading.  These  cracks  will  initiate  at  the  tension  surface  and  continue  upward  to 
approximately  the  neutral  axis.  Analysis  based  upon  this  idea  requires  a different  trans- 
formed cross  section.  This  transformed  cross  section  is  similar  to  that  of  the  uncracked 
transformed  cross  section  except  that  the  concrete  below  the  neutral  axis  is  disregarded. 

Determination  of  the  neutral  axis  for  the  cracked  cross  section  is  similar  to  the  un- 
cracked case.  The  difference  being  that  one  assumes  that  the  concrete  in  tension  does 
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p 


Figure  23:  Compressive  stress  fc  and  tensile  stress  ft  as  a function  of  the  reinforcement 
ratio  p,  based  on  the  uncracked  analysis. 


not  contribute  to  the  moment  calculation.  This  cracked  transformed  section  is  shown 
schematically  in  Fig.  21(c).  The  portion  of  the  concrete  below  the  neutral  axis  is  disre- 
garded, and  the  steel  is  transformed  into  an  equivalent  area  nAs  of  concrete. 

Let  the  depth  of  the  neutral  axis  be  a fraction  A of  the  reinforcement  depth  d.  The 
balance  of  the  area  moments  about  the  neutral  axis  leads  to  the  following  algebraic 
equation: 

b ■■  = nAs(d  — A d)  (88) 

The  solution  for  A is  found  from  the  quadratic  equation: 

A = \J (rip)2  + 2 np  - np  (89) 


The  neutral  axis  depth  factor  A is  shown  in  Fig.  24(a)  as  a function  of  the  reinforce- 
ment ratio  p.  In  contrast  to  the  uncracked  section  analysis,  the  depth  of  the  neutral  axis 
for  the  cracked  cross  section  has  a strong  dependence  on  the  reinforcement  ratio. 

The  stresses  at  the  compression  surface  and  the  steel  reinforcement  are  calculated  in 
much  the  same  way  as  for  the  uncracked  analysis.  The  calculation  of  the  cracked  section 
moment  of  inertia  Icr  omits  the  second  term  on  the  right  hand  side  of  Eqn.  86  because 
this  concrete  is  in  tension.  Assuming  linear  elastic  behavior,  the  compressive  stress  fc  in 
concrete  at  the  top  surface  and  the  tensile  stress  ft  in  the  steel  are  proportional  to  their 
distance  from  the  neutral  axis: 


.fc  = M- 


dX 


.ft  = M. 


d{  1 - A) 


(90) 
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Figure  24:  Cracked  transformed  section  analysis:  (a)  relative  depth  of  neutral  axis  A;  (b) 
relative  tensile  stress  ft/ fy  and  compressive  stress  fc/f'c  for  the  cracked  section  analysis, 
assuming  linear  elastic  behavior. 


The  compressive  stress  fc  and  the  tensile  stress  ft  are  shown  in  Fig.  24(b)  as  a 
function  of  the  reinforcement  ratio  p.  Further,  both  stresses  are  normalized  by  their 
limiting  values.  The  ultimate  compressive  stress  is  ff  and  the  limiting  tensile  stress  is 
the  steel  yield  stress  fy. 

Based  on  the  assumptions  that  the  vault  contents  will  subside  and  that  the  roof  slab 
will  be  a span  in  flexure,  the  top  surface  of  the  slab  should  remain  in  flexure,  regardless 
of  the  shrinkage  strain.  The  effect  of  the  shrinkage  strain,  if  any,  will  be  to  increase  the 
bending  moment  of  the  slab.  Regardless,  the  top  surface  of  a linear  elastic  member  in 
flexure  should  be  in  compression,  with  no  open  cracks.  Therefore,  the  importance  of  the 
structural  analysis  is  to  determine  the  depth  of  the  neutral  axis  and  the  width  of  the 
cracks  within  the  bottom  surface  of  the  slab. 

Proceeding  further  with  the  analysis  requires  employing  a specific  model  for  combined 
cracking  under  flexural  and  shrinkage  stresses.  First,  one  considers  the  development  of 
he  cracks  to  be  independent  of  one  another.  Based  on  this  estimation,  an  estimate  of  the 
distances  between  the  final  cracks  will  be  formulated  based  on  physical  intuition. 

An  estimate  for  the  final  crack  width  and  spacing  within  the  bottom  surface  of  the 
slab  depends  on  the  relative  size  of  the  shrinkage  and  flexural  cracks.  The  estimate  for  the 
ratio  of  shrinkage  crack  spacing  Ls  to  the  flexural  crack  spacing  Lj  is  shown  in  Fig.  25(a) 
as  a function  of  the  reinforcement  ratio  p.  Over  most  of  the  allowable  range  for  p,  the 
ratio  Ls/Lf  is  greater  than  five.  Therefore,  the  flexural  cracks  should  develop  between 
shrinkage  cracks. 

An  independent  analysis  of  the  shrinkage  crack  width  ws  and  flexural  crack  width  Wf 
are  shown  in  Fig.  25(b)  as  a function  of  the  reinforcement  ratio  p.  From  an  assessment 
of  the  estimate  of  crack  widths,  the  combined  crack  widths  will  certainly  be  tens,  if  not 
hundreds,  of  micrometers  wide.  The  cracks  will  occur  at  the  tension  face  and  extend  at 
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Figure  25:  Results  of  the  cracking  model:  (a)  The  ratio  of  the  shrinkage  crack  spacing 
Ls  and  the  flexural  crack  spacing  Lf  \ (b)  Flexural  crack  width  Wf  and  shrinkage  crack 
width  ws. 


least  to  the  neutral  axis. 


7.4.9  Simplified  Composite  Permeability 

Since  the  permeability  of  the  cracked  portion  of  the  slab  will  be  orders  of  magnitude 
greater  than  the  uncracked  portion  of  the  concrete,  further  analysis  can  be  simplified  by 
approximating  the  composite  permeability  k*  using  Eqn.  56: 


This  result  is  plotted  in  Fig.  26(a)  as  a function  of  the  reinforcement  ratio  p.  The  ratio  of 
the  permeabilities  increases  by  less  than  a factor  of  five  over  most  of  the  allowable  range 
of  p. 

Since  the  exact  crack  width  at  the  bottom  surface  of  the  concrete  may  become  in- 
consequential, the  significance  of  the  structural  analysis  will  be  to  formulate  an  accurate 
estimate  for  the  depth  of  the  neutral  axis.  This  fact  is  emphasized  in  Fig.  26(b)  that 
shows  a plot  of  the  ratio  k* /kc  as  a function  of  the  relative  depth  to  the  neutral  axis  a,  for 
three  values  of  the  ratio  kb/kc.  Note  that  the  ratio  k* /kc  converges  rapidly  with  respect 
to  the  ratio  kb/kc.  In  practice,  the  ratio  kb/kc  could  be  orders  of  magnitude  greater  than 
those  shown  in  the  graph. 


7.5  Probabilistic  Approach 

A deterministic  calculation,  like  that  performed  in  the  previous  section,  is  useful  for  esti- 
mating the  probable  behavior  of  the  reinforced  concrete  structure.  In  practice,  however, 
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Figure  26:  Ratio  of  the  composite  permeability  k*  to  the  uncracked  concrete  permeability 
kc:  (a)  as  a function  of  reinforcement  ratio  p,  and  (b)  as  a function  of  the  relative  depth 
to  the  neutral  axis  a. 

the  exact  size  of  the  cracks  that  form  will  vary,  as  will  the  depth  of  the  neutral  axis.  The 
resulting  variations  in  these  quantities  lead  to  variations  in  the  expected  permeability  of 
the  structure.  In  order  to  capture  this  variation  a probabilistic  approach  is  used. 

The  probabilistic  method  employed  here  is  based  on  a Monte  Carlo  calculation  using 
parameter  uncertainty  as  a means  of  characterizing  random  deviates.  Each  parameter 
required  for  a deterministic  calculation  has  associated  with  it  an  estimated  uncertainty. 
With  this  uncertainty,  one  can  use  a computer  to  repeat  the  deterministic  calculation  nu- 
merous times,  with  each  iteration  using  a new  random  deviate  for  each  parameter.  From 
the  multiple  calculations,  the  estimated  composite  permeability  from  the  each  iteration 
form  a population  of  permeability  values,  from  which  a population  mean  and  standard 
deviation  can  be  calculated.  These  values  can  then  be  incorporated  into  the  service  life 
calculation,  which  will,  in  turn,  yield  an  estimated  mean  and  standard  deviation  for  the 
service  life  of  the  structure. 

7.5.1  Initial  Material  Parameters 

The  following  probabilistic  analysis  follows  identically  from  the  previous  section  with  the 
only  exception  being  that  some  of  the  input  parameters  are  normally  distributed  about 
some  mean  value  with  a coefficient  of  variation  (CoV)  (the  ratio  of  the  mean  to  the 
standard  deviation).  The  parameters  allowed  to  vary  are  shown  in  Table  10.  Note  that 
the  mean  values  are  identical  to  those  used  in  the  previous  deterministic  analysis.  The 
values  for  the  (CoV)  are  not  meant  to  be  accurate  for  any  particular  structure,  and  are 
only  given  for  demonstration  of  the  technique.  The  reader  is  referred  to  the  ACI  117 
report  on  construction  tolerances  for  general  guidance. 

The  values  for  the  coefficient  of  variation  given  in  Table  10  are  not  based  on  any 
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Table  10:  The  mean  and  associated  coefficient  of  variation  (CoV)  for  the  input  parameters 
that  would  be  used  in  an  uncertainty  analysis  of  flexural  and  shrinkage  cracking. 


Variable 

Property 

Mean 

CoV 

L 

Slab  length  (m) 

10.0 

0.05 

Ec 

Concrete  Modulus  (GPa) 

28.6 

0.05 

f'c 

Concrete  compressive  strength  (MPa) 

35.0 

0.05 

h 

Slab  thickness  (m) 

1.00 

0.10 

d 

Depth  of  steel  (in) 

0.90 

0.05 

Shrinkage  strain  (pm/m) 

90 

0.25 

particular  data.  An  analysis  of  particular  vault  should  be  done  with  careful  consideration 
of  the  construction  plans  and  the  project  description.  For  example,  the  CoV  for  the  slab 
thickness  given  above  is  large  compared  with  normal  slab  on  grade  construction.  The 
reason  for  this  is  the  assumption  that  the  contents  of  the  vault  would  be  compacted 
remotely,  arid  so  the  variability  of  the  slab  thickness  would  be  greater  than  for  normal 
practices.  Construction  plans  for  specific  vaults  may  vary  and  warrant  a smaller  CoV. 

For  the  probabilistic  analysis,  the  randomization  is  done  as  a function  of  the  designed 
reinforcement  ratio  p.  For  each  value  of  p,  random  values  are  drawn  for  the  parameters 
shown  above.  Using  these,  the  analysis  proceeds  as  for  the  deterministic  analysis.  The 
approach  used  by  4SIGHT  would  be  to  establish  a fixed  reinforcement  ratio  p and  then 
repeat  the  analysis  at  a constant  value  of  p,  randomizing  the  input  parameters  listed 
above.  4SIGHT  would  then  report  the  mean  and  standard  deviation  of  the  results  at  a 
given  reinforcement  ratio. 

7,5.2  Results 

The  material  and  structural  parameters  given  in  Table  10  are  used  to  recalculate  the 
same  quantities  as  were  calculated  for  the  deterministic  analysis.  The  only  difference  is 
that  the  calculation  is  repeated  many  times.  Each  time,  a value  for  the  reinforcement 
ratio  p is  chosen  and  then  the  computer  generates  a random  number  for  each  parameter 
listed  in  Table  10.  These  “randomized"  values  are  then  used  in  the  deterministic  analysis. 
The  variation  in  the  result  indicates  the  output  sensitivity  to  that  parameter,  at  that 
value  of  CoV.  Each  figure  of  the  results  to  be  shown  includes  a solid  curve  indicating  the 
values  from  the  previous  deterministic  calculations. 

The  results  for  the  depth  of  the  neutral  axis  and  the  relative  stresses  are  shown  in 
Fig.  27  and  can  be  compared  directly  to  Fig.  24.  In  general,  the  probabilistic  data  follow 
the  same  trends  as  in  the  deterministic  case.  The  calculated  depths  of  the  neutral  axis 
shown  in  Fig.  27(a)  are  a function  of  the  reinforcement  ratio  p.  The  results  indicate  that 
the  uncertainty  in  the  depth  of  the  neutral  axis  is  relatively  insensitive  to  the  uncertainty 
in  the  input  parameters. 

By  comparison,  the  relative  compressive  and  tensile  stresses  shown  in  Fig.  27(b)  are 
much  more  sensitive.  The  variability  in  the  compressive  stress  ratio  fc/ f'c  is  dramatic, 
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Figure  27:  Probabilistic  analysis  results:  (a)  Relative  depth  of  neutral  axis  A;  (b)  Relative 
tensile  stress  ft/fy  and  compressive  stress  fc/ f'c  for  the  cracked  section  analysis. 


with  some  of  the  values  approaching  one,  even  for  relatively  large  reinforcement  ratios. 
Since  the  analysis  has  been  done  assuming  a linear  elastic  material,  a relative  stress 
greater  than  0.5  is  probably  unreliable.  Regardless,  the  tensile  stresses  in  the  bottom  of 
the  slab  will  control  the  crack  width. 

The  ratio  of  the  shrinkage  crack  spacing  Ls  to  the  flexural  crack  spacing  Lf  is  shown 
in  Fig.  28(a)  as  a function  of  the  reinforcement  ratio  p.  Although  the  ratio  has  large 
variation,  the  value  is  typically  far  greater  than  one.  Therefore,  regardless  of  the  variation 
in  input  parameters,  one  expects  that  flexural  cracks  will  be  initiated  between  the  existing 
shrinkage  cracks. 

The  shrinkage  crack  width  ws  and  flexural  crack  width  wj  are  shown  in  Fig.  28(b).  The 
shrinkage  crack  widths  seem  to  be  insensitive  to  the  variations  in  the  input  parameters. 
Although  the  flexural  cracks  appear  to  be  more  sensitive,  the  data  are  plotted  on  a 
logarithmic  ordinate  which  enhances  the  appearance  of  variability  in  smaller  quantities. 

Since  the  expected  crack  widths  are  still  on  the  order  of  tens  or  hundreds  of  microm- 
eters, the  estimate  of  the  bulk  concrete  permeability  k*  can  again  be  simplified  with  the 
following  approximation: 

k*  h 
kc  A d 

This  result  is  shown  in  Fig.  29.  Since  the  depth  of  the  neutral  axis  was  fairly  insensitive 
to  the  variability  of  the  input  parameters,  the  relative  bulk  permeability  is  also  relatively 
insensitive  to  parameter  variation.  The  added  variation  over  A is  due  to  the  variability 
in  the  roof  thickness  h and  the  depth  of  steel  reinforcement  d. 

The  result  shown  in  Fig.  29  is  significant  because  it  reveals  that  the  details  of  the 
drying  shrinkage  model  are  relatively  unimportant.  Rather,  the  most  significant  factors 
in  predicting  the  permeability  of  a reinforced  concrete  structure  is  the  permeability  of  the 
sound  concrete,  and  the  depth  of  the  neutral  axis.  Therefore,  the  two  most  important 
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Figure  28:  Results  of  the  probabilistic  analysis  of  the  cracking  models:  (a)  the  ratio  of 
the  shrinkage  crack  spacing  Ls  to  the  flexural  crack  spacing  Lf  \ (b)  flexural  crack  width 
Wf  and  shrinkage  crack  width  ws. 


requirements  of  the  engineering  using  4SIGHT  to  assess  a concrete  vault  is  an  accurate 
estimation  of  the  sound  concrete  permeability  and  an  accurate  estimate  of  the  depth  of 
the  neutral  axis. 
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Figure  29:  From  the  probabilistic  analysis,  the  ratio  of  the  composite  permeability  k*  to 
the  uncracked  concrete  permeability  kc. 
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8 User’s  Guide 


An  evaluation  of  the  previous  user  interface  to  4SIGHT  revealed  unnecessary  complex- 
ity. In  an  effort  to  simplify  the  use  of  4SIGHT,  a number  of  assumptions  have  been 
incorporated.  The  input  file  is  now  more  restrictive,  but  the  requirements  are  fewer.  The 
objective  was  to  simplify  the  user’s  burden  and  facilitate  the  independent  development  of 
a graphical  front  end  for  user  input  and  computed  output.  A graphical  front  end  could 
be  tailored  to  a specific  software  platforms,  while  the  computational  portion  remains 
relatively  platform  independent. 

8.1  Input  File 

In  general,  most  of  the  parameters  from  4SIGHT  remain.  In  some  cases,  however,  the 
expression  has  been  compressed  to  facilitate  the  addition  of  parameter  uncertainty.  A 
number  of  parameters  require  a single  value,  others  require  the  specification  of  a proba- 
bility distribution  (for  parameter  uncertainty  Monte  Carlo  calculation),  and  still  others 
allow  for  default  values  with  a specified  coefficient  of  variation.  A useful  means  of  de- 
scribing the  input  file  is  to  give  an  example: 

21  N_cells 
4 num_species 
8 0.000  0.300  0.000 

19  0.000  0.200  0.000 

11  0.200  0.100  0.000 
17  0.200  0.000  0.000 

0 1 f ree_surf ace_bool 
0.02  0.0  pressure_MPa 
-4  -9.00  0.500  permeability_cm/s 
-3  1000.  0.25  formation.f actor 
-3  0.150  0.05  porosity 
+3  0.500  0.05  wc 
+3  20.0  1.00  f c_MPa 
-3  20.0  0.10  Elastic_modulus_GPa 
3 10.0  0.10  length_m 

3 1.000  0.020  thickness.m 

3 0.900  0.010  steel_depth_m 

3 10.00  1.0  soil_depth_m 

-3  0.010  0.050  As/Ac  : remf  orcement_rat  10 
100000.  time_limit_days 
23204  seed 
10  simulations 

The  parameters  shown  will  be  discussed  on  a line  by  line  basis.  The  optional  text  that 
follows  each  numerical  parameter(s)  is  included  as  an  aid  to  the  user,  but  is  ignored  by 
the  computer  program. 
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Line  1:  N_cells 

The  number  of  finite  difference  computational  elements  used  to  represent  the  slab  for  the 
1-D  calculation.  A value  of  21  is  typically  sufficient. 

Line  2:  num.species 

The  number  of  ionic  species  available  for  transport.  The  available  ionic  species  are  the 
following: 


Key 

Species 

1 

H+ 

3 

Li+ 

11 

Na+ 

12 

Mg+2 

19 

K+ 

20 

Ca+2 

6 

C03~2 

7 

N03- 

8 

OH- 

16 

S04-2 

17 

CD 

35 

Br- 

53 

I- 

In  this  example  there  are  four  species:  OH  , Cl  , K+,  Na+. 

Line  2a, 2b,. . . : 

The  quantity  of  each  of  the  mim_species  present.  The  four  values  given  (AtNo  cA  clnt 
cB)are  as  follows: 

AtNo  : atomic  number  of  species.  In  the  case  of  OH-  and  SO;j-,  it  is  the  atomic  number 
of  the  first  species.  Ex.:  OH-  = 8,  SO4-  = 16,  etc. 

cA  : molarity  of  species  external  to  Side  A. 

clnt  : molarity  of  species  within  the  concrete. 

cB  : molarity  of  species  external  to  Side  B. 

Side  A and  Side  B are  arbitrary. 

Line  3:  f ree_surf ace_bool 

A boolean  value,  either  0 or  1,  denoting  whether  Side  A or  Side  B is  a f ree_surf  ace. 
For  a f ree_surf  ace,  the  diffusive  flux  is  zero  and  ions  within  the  pore  solution  accu- 
mulate m the  presence  of  a Darcy  flow. 

Line  4:  pressure_MPa 

The  hydrostatic  pressure  at  side  0 and  side  i,  in  units  of  MPa. 
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n.b.:  For  some  of  the  following  parameters,  the  user  is  to  supply  4SIGHT  with  information 
regarding  the  type  of  distribution  to  use  and  the  parameters  of  that  distribution:  xQ  and 
a.  The  meaning  of  the  distribution  parameters  have  been  discussed  previously.  The 
particular  distribution  is  specified  using  the  following  numbering  scheme: 

n = 1 : delta  function 

n = 2 : uniform  distribution 

n = 3 : Gaussian  distribution 

n = 4 : Lognormal  distribution 

If  the  user  is  specifying  a definite  value  (and  not  reiving  upon  4SIGHT  to  create  and 
estimated  value  for  them),  the  distribution  type  n is  made  positive  using  a leading  plus 
sign.  If  the  distribution  type  n is  made  negative  with  a leading  minus  sign,  4SIGHT 
will  generate  a default  value  based  on  the  remaining  specified  values.  If  a default  value 
is  requested,  the  value  for  xQ  appearing  in  the  input  file  will  be  ignored  and  and  the 
user  specifies  the  coefficient  of  variation  of  this  default  value  through  the  parameter 
a.  The  exception  to  this  is  the  hydraulic  conductivity:  if  its  distribution  type  n is 
negative,  o represents  the  standard  deviation  of  the  logarithm,  which  is  identical  to  its 
meaning  for  positive  values  of  the  distribution  type.  The  relationships  used  by  4SIGHT 
to  calculate  default  values  either  appear  in  this  report  or  can  be  found  in  the  previous 
4SIGHT  report  [1], 

Line  5:  permeability_cm/s 

The  logarithm  (base  10)  of  the  hydraulic  conductivity  of  uncracked  concrete,  when  ex- 
pressed in  units  of  cm/s.  The  three  values  shown  are  the  distribution  type,  xQ,  and  a. 
In  the  example,  the  user  is  requesting  that  the  computer  estimate  the  hydraulic  conduc- 
tivity, and  that  the  standard  deviation  of  the  logarithms  (base  10)  of  these  values  be 
0.500.  This  corresponds  to  hydraulic  conductivity  values  varying  over  one  half  and  order 
of  magnitude. 

Line  6:  formation_f actor 

The  formation  factor  of  the  intact  concrete.  The  three  numbers  are  the  distribution  type, 
xQ,  and  a.  In  the  example  above,  the  user  is  requesting  4SIGHT  to  estimate  the  formation 
factor  from  the  other  parameters,  and  that  these  predicted  values  will  have  a coefficient 
of  variation  of  25  %. 

Line  7:  porosity 

The  capillary  porosity  of  the  intact  concrete.  The  three  numbers  are  the  distribution 
type,  r0,  and  a.  In  the  example  above,  the  user  is  requesting  4SIGHT  to  estimate  the 
porosity  from  the  other  parameters,  and  that  these  estimated  values  will  have  a coefficient 
of  variation  of  5 %. 
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Line  8:  wc 

The  water  cement  ratio  at  the  time  of  mixing.  The  three  numbers  are  the  distribution 
type,  xD1  and  a.  In  the  example  above,  the  user  is  specifying  that  the  water  cement  ratio 
shall  be  chosen  from  a Gaussian  distribution  with  mean  0.50  and  standard  deviation  0.05. 

Line  9:  fc_MPa 

The  design  compressive  strength  of  the  concrete.  The  three  numbers  are  the  distribution 
type,  xa,  and  a.  In  the  example  above,  the  user  is  specifying  that  the  compressive 
strength  shall  be  chosen  from  a Gaussian  distribution  with  mean  20  MPa.  and  standard 
deviation  1.0  MPa. 

n.b.:  At  least  one  of  the  above  six  parameters  must  be  specified  by  the  user.  Therefore, 
at  least  one  of  these  values  must  have  a positive  distribution  type.  Otherwise,  4SIGHT 
has  no  information  from  which  to  characterize  the  concrete.  Of  course,  the  user  may 
specify  all  of  these  parameters.  In  the  above  input  file,  the  + and  - are  used  to  remind 
the  user  of  this  requirement. 

Line  10:  Elastic_modulus_GPa 

The  elastic  modulus  of  the  concrete,  in  units  of  GPa.  The  three  numbers  are  the  distribu- 
tion type,  xQ1  and  a . In  the  example  above,  the  user  is  specifying  that  the  elastic  modulus 
be  estimated  from  the  other  material  parameters,  and  that  these  values  be  chosen  to  have 
a coefficient  of  variation  of  10  %. 

Line  11:  length_m 

The  length  of  the  slab,  in  units  of  m.  The  three  numbers  are  the  distribution  type,  xQl 
and  a.  In  the  example  above,  the  user  is  specifying  that  the  length  be  estimated  from  a 
Gaussian  distribution  with  mean  10.00  and  standard  deviation  0.10. 

Line  12:  thickness.m 

The  vertical  thickness  of  the  slab,  in  units  of  m.  The  three  numbers  are  the  distribution 
type,  xQ,  and  a.  In  the  example  above,  the  user  is  specifying  that  the  thickness  be 
estimated  from  a Gaussian  distribution  with  mean  1.00  and  standard  deviation  0.02. 

Line  12:  steel_depth_m 

The  depth  of  reinforcement  steel,  in  units  of  m.  The  three  numbers  are  the  distribution 
type,  xQ , and  a.  In  the  example  above,  the  user  is  specifying  that  the  depth  of  steel 
reinforcement  be  estimated  from  a Gaussian  distribution  with  mean  0.90  and  standard 
deviation  0.01. 

Line  12:  soil_depth_m 

The  depth  of  soil  over  the  vault,  in  units  of  m.  The  three  numbers  are  the  distribution 
type,  x0 , and  a.  In  the  example  above,  the  user  is  specifying  that  the  depth  of  steel 
reinforcement  be  estimated  from  a Gaussian  distribution  with  mean  10.00  and  standard 
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deviation  1.00. 

Assigning  a negative  value  to  the  depth  of  soil  informs  4SIGHT  that  the  user  does  not 
want  a crack  analysis,  and  that  the  concrete  should  remain  entirely  intact,  without  any 
cracks. 

Line  12:  As/Ac  : remf  orcement_ratio 

The  reinforcement  ratio,  defined  as  the  ratio  of  the  area  of  steel  (As)  to  the  area  of  concrete 
(Ac).  The  three  numbers  are  the  distribution  type,  xQ , and  a.  In  the  example  above, 
the  user  is  specifying  that  4SIGHT  should  estimate  a “best”  guess  from  the  available  data 
and  that  deviates  should  have  a coefficient  of  variation  of  5 %. 

Line  12:  time_limit_days 

The  time  limit  for  the  service  life  calculation,  in  units  of  days. 

Line  12:  seed 

A seed  for  the  random  number  generator.  Using  the  same  seed  will  result  in  identical 
values  being  generated. 

Line  12:  simulations 

The  number  of  times  to  perform  a parameter  uncertainty  calculation. 
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9 Summary 

These  efforts  directed  towards  4SIGHT  have  moved  it  towards  a rigorous  model  of  advective- 
diffusive  transport  within  concrete.  This  is  a fundamental  part  of  establishing  a service 
life  model  based  on  the  simulation  of  reactions  within  the  concrete.  Every  degradation 
mechanism  relevant  to  concrete  performance  depends,  in  some  way,  on  the  transport 
characteristics  of  the  material.  Although  applying  transport  models  to  both  reference 
data  and  carefully  controlled  experiments  appear  to  be  logical  endeavors,  they  are  rare 
procedures  in  the  concrete  materials  research  literature. 

The  work  done  for  this  project  has  shown  that  a carefully  constructed  model  that 
accounts  for  ion-ion  interactions,  due  to  both  electrostatic  interactions  and  excluded 
volume  interactions,  can  then  be  used  to  model  “real  world”  behavior.  Further,  the 
resulting  transport  equation  reveals  the  meaningful  physical  quantities  that  are  required 
from  experiments  performed  on  the  material  of  interest.  It  was  serendipitous  that  the 
more  rigorous  equation  actually  requires  fewer  physical  parameters,  regardless  of  the 
number  of  different  species  being  transported. 

An  unforeseen  result  of  this  work  has  been  a critical  evaluation  of  the  ineaningful- 
ness  of  the  apparent  diffusion  coefficient.  Experiments  verified  theoretical  expectations 
that  the  apparent  diffusion  coefficient  can  vary  over  time  in  nonreacting  porous  systems. 
Data  from  experiments  done  here  have  shown  that  the  apparent  diffusion  coefficient  of  a 
particular  ion  can  vary  by  a factor  of  ten,  depending  on  the  other  ionic  species  present. 
Further,  a particular  collection  of  species  can  lead  to  a negative  apparent  diffusion  coef- 
ficient. In  contrast,  a rigorous  transport  equation,  using  a fixed  porosity  and  formation 
factor,  can  capture  this  behavior  by  incorporating  the  ion-ion  interactions  into  the  trans- 
port equation. 

The  addition  of  a simple  structural  analysis  into  4SIGHT  will  allow  users  to  make,  in 
the  absence  of  a more  thorough  analysis,  reasonable  approximations  of  the  width,  depth, 
and  spacing  of  drying  shrinkage  and  flexural  cracks.  An  analytical  investigation  of  these 
relationships  revealed  that  the  most  important  factors  influencing  the  bulk  permeability 
of  a reinforced  concrete  slab  in  flexure  is  the  permeability  of  the  uncracked  concrete 
and  the  depth  of  the  neutral  axis.  The  permeability  is  insensitive  to  the  details  of  the 
particular  crack  model  employed. 

The  addition  of  a Monte  Carlo  calculation  based  on  parameter  uncertainty  addresses 
the  inherent  issue  of  variability.  A deterministic  service  life  calculation  is  meaningless 
since  it  would  be  impossible  to  guarantee  the  predicted  result.  Rather,  by  at  least 
incorporating  the  parameter  uncertainty  into  the  result,  one  can  make  a statement  about 
the  result  uncertainty  as  a function  of  the  parameter  uncertainty.  Further,  a probabilistic- 
result  is  a better  representation  of  a physical  problem  since  one  can  then  ask  for  prediction 
intervals  of  the  computed  results. 
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